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We identify 31 dimensionless physical constants required by particle physics and cosmology, and 
emphasize that both microphysical constraints and selection effects might help elucidate their origin. 
Axion cosmology provides an instructive example, in which these two kinds of arguments must both 
be taken into account, and work well together. If a Peccei-Quinn phase transition occurred before or 
during inflation, then the axion dark matter density will vary from place to place with a probability 
distribution. By calculating the net dark matter halo formation rate as a function of all four relevant 
cosmological parameters and assessing other constraints, we find that this probability distribution, 
computed at stable solar systems, is arguably peaked near the observed dark matter density. If 
cosmologically relevant WIMP dark matter is discovered, then one naturally expects comparable 
densities of WIMPs and axions, making it important to follow up with precision measurements to 
determine whether WIMPs account for all of the dark matter or merely part of it. 



I. INTRODUCTION 

Although the standard models of particle physics and 
cosmology have proven spectacularly successful, they to- 
gether require 31 free parameters (Table 1). Why we 
observe them to have these particular values is an out- 
standing question in physics. 

A. Dimensionless numbers in physics 

This parameter problem can be viewed as the logical 
continuation of the age-old reductionist quest for simplic- 
ity. Realization that the material world of chemistry and 
biology is built up from a modest number of elements 
entailed a dramatic simplification. But the observation 
of nearly one hundred chemical elements, more isotopes, 
and countless excited states eroded this simplicity. 

The modern 5C/(3) x SU{2) x U{1) standard model 
of particle physics provides a much more sophisticated 
reduction. Key properties (spin, electroweak and color 
charges) of quarks, leptons and gauge bosons appear as 
labels describing representations of space-time and in- 
ternal symmetry groups. The remaining complexity is 
encoded in 26 dimensionless numbers in the Lagrangian 
(Table 1)"'^. All current cosmological observations can be 
fit with 5 additional parameters, though it is widely antic- 
ipated that up to 6 more may be needed to accommodate 
more refined observations (Table 1). 

Table 2 expresses some common quantities in terms of 
these 31 fundamental ones^, with ~ denoting cruder ap- 
proximations than Many other quantities commonly 



^ Here and A are defined so that the Higgs potential is V("I>) = 

/.2|$|2 + A|$|4. 
^ The last six entries are mere order-of- magnitude estimates 

In the renormalization group approximation for a in Table 2, 



referred to as parameters or constants (see Table 3 for a 
sample) are not stable characterizations of properties of 
the physical world, since they vary markedly with time 
0. For instance, the baryon density parameter fif,, the 
baryon density pb, the Hubble parameter h and the cos- 
mic microwave background temperature T all decrease 
toward zero as the Universe expands and are, de facto, 
alternative time variables. 

Our particular choice of parameters in Table 1 is a 
compromise balancing simplicity of expressing the funda- 
mental laws (i.e., the Lagrangian of the standard model 
and the equations for cosmological evolution) and ease of 
measurement. All parameters except /z^, pA, ^b, and 

are intrinsically dimensionless, and we make these final 
five dimensionless by using Planck units (for alternatives, 
see 0, 0)- Throughout this paper, we use "extended" 
Planck units defined by c = G = ?i = |ge| = fcs = 1- We 
use fi — 1 rather than ft, = 1 to minimize the number of 
(27r)-factors elsewhere. 



B. The origin of the dimensionless numbers 

So why do we observe these 31 parameters to have the 
particular values listed in Table 1? Interest in that ques- 
tion has grown with the gradual realization that some 
of these parameters appear fine-tuned for life, in the 
sense that small relative changes to their values would 
result in dramatic qualitative changes that could pre- 
clude intelligent life, and hence the very possibility of 
refiective observation. As discussed extensively elsewhere 
iS, da, dL 12, ii, I£ IL JJ. 2Q, ,2i, 2aj , there 
are four common responses to this realization: 



only fermions with mass below mz should be included. 
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Table 1: Input physical parameters." Those for particle physics (the first 26) and pA appear explicitly in the the Lagrangian, whereas 
the cosmological ones (the last 11) are inserted as initial conditions. The last six are currently optional, but may become required to fit 
improved measurements. The values are computed from the compilations in |^|^|^. 



Parameter 


Meaning 


Measured value 


g 


Weak coupling constant at mz 


0.6520 ±0.0001 




Weinberg angle 


0.48290 ± 0.00005 


9a 


Strong coupling constant at mz 


1.221 ± 0.022 




Quadratic Higgs coefficient 


10-^^ 


A 


Quartic Higgs coefficient 


~ 1? 


Ge 


Electron Yukawa coupling 


2.94 X 10-^ 


Gy. 


Muon Yukawa coupling 


0.000607 


Gr 


Tauon Yukawa coupling 


0.0102156233 


Gu 


Up quark Yukawa coupling 


0.000016 ± 0.000007 


Gd 


Down quark Yukawa coupling 


0.00003 ± 0.00002 


Gc 


Charm quark Yukawa coupling 


0.0072 ± 0.0006 


Gs 


Strange quark Yukawa coupling 


0.0006 ± 0.0002 


Gt 


Top quark Yukawa coupling 


1.002 ±0.029 


Gb 


Bottom quark Yukawa coupling 


0.026 ± 0.003 


sin 012 


Quark CKM matrix angle 


0.2243 ± 0.0016 


sin 023 


Quark CKM matrix angle 


0.0413 ± 0.0015 


sin 013 


Quark CKM matrix angle 


0.0037 ± 0.0005 


Sl3 


Quark CKM matrix phase 


1.05 ± 0.24 


6qcd 


CP-violating QCD vacuum phase 


< 10-9 


Gve 


Electron neutrino Yukawa coupling 


< 1.7 X 10-" 


Gvfj, 


Muon neutrino Yukawa coupling 


< 1.1 X 10-** 


G.^ 


Tau neutrino Yukawa coupling 


< 0.10 


sin 0'-^2 


Neutrino MNS matrix angle 


0.55 ± 0.06 


sin 2^23 


Neutrino MNS matrix angle 


> 0.94 


sin 0'-^^ 


Neutrino MNS matrix angle 


< 0.22 


^'la 


Neutrino MNS matrix phase 


7 


PA 


Dark energy density 


(1.25 ±0.25) X 10-123 


€b 


Baryon mass per photon ph/i^-y 


(0.50 ±0.03) X 10~28 




Cold dark matter mass per photon pc/nj 


(2.5 ± 0.2) X 10-2** 




Neutrino mass per photon pvjn^ = X] 


< 0.9 X 10-28 


Q 


Scalar fluctuation amplitude <5h on horizon 


(2.0 ± 0.2) X 10-5 


ris 


Scalar spectral index 


0.98 ± 0.02 


otn 


Running of spectral index dris/dlnfc 


|a| < 0.01 


r 


Tensor-to-scalar ratio (Qt/Q)^ 


<0.36 


nt 


Tensor spectral index 


Unconstrained 


w 


Dark energy equation of state 


-1 ± 0.1 


K 


Dimensionless spatial curvature k/a?T'^ [85] 


\k\ < 10-60 



"Throughout this paper, we use extended Planck units c = h = 
G = fcf, = |5e| = 1. For reference, this common convention gives 
the following units: 



Constant 


Definition 


Value 


Planck length 
Planck time 
Planck mass 
Planck temperature 
Planck energy 
Planck density 
Unit charge 
Unit voltage 


(nG/c3)l/2 
(?iG/c5)l/2 

(ac/G)i/2 
{hey ay ^^/k 

c^/hG'^ 

\le\ 

(?icVG)l/2/|g,| 


1.61605 X lO-^'^m 
5.39056 X 10-*-*s 
2.17671 X 10-*kg 
1.41696 X 10^2^ 
1.22105 X lO^^GeV 
5.15749 X lO^'^kg/m^ 
1.60218 X 10-19 C 
1.22105 X 102* V 
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lable i: Derived physical paramet 


ers, 


in extended Planck units c = h = G = k), 


1 1 -1 
= kel = 1. 




Cleaning 


Definition 


IVIeasurecl value 


e 


Electromagnetic coupling constant at mz 




Ri 0.313429 ± 0.000022 




Electromagnetic interaction strength at mz 


e2/47r = sin^ 6'w/47r 


K. 1/(127.918 ± 0.018) 


Ctw 


Weak interaction strength at mz 




~ 0.03383 ± 0.00001 


Cla 


Strong interaction strength at mz 




~ O.II8D ± 0.0042 


ag 


Gravitational couphng constant 


Gm^/hc = rrip 


~ 5.9046 X 10 

-1 


a 


Electromagnetic interaction strength at 




1 \ 1 , 2 1 m^^ 


1/137.03599911(46) 




aym-A ^ + 7^- In — j — , 3 — , — 3- 




mass 


vg/2 


(80 425 + 0381 GeV 




Z mass 


vg/2 cos By^r 


(91.1876 ± 0.0021) GeV 


"J* 


Fermi constant 




Si 1.17 X 10~5 GeV~2 


'"'±1 


Higgs mass 


V- 


-p2/2 


100-250 GeV? 


V 


Higgs vacuum expectation value 


V- 




(246.7 ± 0.2) GeV 


trie 


Electron mass 


vGe/V2 


(510998.92 ± 0.04)eV 




]VIuon mass 


vGf,/V2 


(105658369 ± 9)eV 




Tauon mass 


vG-r/V2 


(1776.99 ± 0.29) MeV 


TJT/y, 


Up c[uark mass 


vGu/V2 


(1.5 — 4) MeV 


•"'a 


Down c[uark mass 


vGg/V2 


(4 - 8) MeV 


TTT-c 


Charm Quark mass 


vGc/V2 


CI 15 — 1 351 GeV 




Strange cj^uark mass 


vGs/V2 


I'sn _ isnl MeV 


TtT/f 


' 1 V\Y^ nAM^vVc rVlflQQ 
_L W LJ ^ LiCLi JV lllOjOO 


vGt/V2 


(■1 74 Q -L K 1 ^ OeV 


TYl}f 


Bottom quark mass 


vGb/V2 


(4,1 _ 4.9) GeV 






vG^jV2 


< 3eV 


Tn.,, 


JVIuon neutrino mass 


vG,^/V2 


< 0.19 MeV 




-LOiLi. llCLiUllllU Hldoo 


vG^jV2 


< 1 8 2 GeV 


TTlp 


Proton mass 


2m 


„ + md+QCD+QED 


(938.27203 ± 0.00008) MeV 


mn 


Neutron mass 


2ma + mu+QCD+QED 


(939.56536 ± 0.00008) MeV 




Electron/proton mass ratio 


me/mp 


1/1836.15 


I3„ 


Neutron/proton relative mass difference 


m„ 


Imp - 1 


1/725.53 


Ry 


Hydrogen binding energy (Rydberg) 


Ry 


= m^c^c? 12 = o2/3mp/2 


K. 13.6057eV 


Ob 


Bohr radius 


as 


= h/cmea = {al3mp)~^ 


ss 5.29177 X 10-"m 


o-t 


Thomson cross section 


Sir 
3 




Ri 6.65246 X 10-29 


kc 


Coulomb's constant 


l/47reo = hca/q'^ = a 


1/137.03599911(46) 


V 


Baryon/photon ratio 


rih/rij =^i,/mp 


(6.3 ±0.3) X 10-1° 




Matter per photon 




f ?c + = pm/n^ = m^Tjil + Rc + Rv) 


(3.3 ±0.3) X 10-28 ~ 4cV 




CDM/baryon density ratio 


Pc/Pb = ?c/?b = '^cdm/'^b 


K. 6 


Ri/ 


Neutrino/baryon density ratio 


/ fv^ai 3 
PvlPh — 11 ^ 


rZ 1 


Mu 


Sum of neutrino masses 




= Zpulriv = {ll/3)mp77_Ri/ 






Neutrino density fraction 




< 0.1 




Matter-radiation equality temperature 


30C(3)N + 21 (4)4/31 i^^o.220189C 


Ki 9.4 X lO^K 


eq 


Matter density at equality 


Pm 


;Teq) = 765314352000C(3)'^ ,4 ^ 0.00260042^* ^ 1.1 X lO-i^kg/m^ 

^ ^' (242+21x222/3)^7rl<l^ ^ 




Dark energy domination epoch 


= (PmVPA)i/« ~ 0.137514e4/3p;^i/3 


3215 ± 639 


'"'galaxy 


Mass of galaxy ^] 


~ a5^-l/2^-3 


~ 10"! kg 


TYlit 


Mass of star [21 

mm 


-2 
~ rrip 


~ lO^Okg 


'^planet 


Mass of habitable planet [?| 




~ 1024kg 


'^asteroid 


Maximum mass of asteroid [2] 


~ a3/2^3/2„-i/2 


~ 10^2 kg 


'^person 


Maximum mass of person [2] 


~ 102a3/4^-l/2 


~ 102kg 


^wimp 


WIMP dark matter density per photon 


~ la^v'^/g'^ = -lOVW 


~ 3 X 10-28? 



1. Fluke: Any apparent fine-tuning is a fluke and is 
best ignored. 

2. Multiverse: These parameters vary across an en- 
semble of physically realized and (for all practical 



purposes) parallel universes, and we flnd ourselves 
in one where life is possible. 

3. Design: Our universe is somehow created or sim- 
ulated with parameters chosen to allow life. 
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Table 3: Derived physical variables, in extended Planck units c = h = G = ki, = \qe\ = I. 



Parameter 


Meaning 


Definition 


Measured value 


T 


CMB temperature 


(Acts as a time variable) 


(2.726 ± 0.005)K (today) 


n^y 


Photon number density 


^^T^, C(3) =i 1.20206 


0.243588r3 


riv 


Neutrino number density 


9„ _ 18C(3)^3 
11 ' IItt^ 


0.199299T3 


r y 


Photon density 


15 


0.657974r'' 


Pb 


Baryon density 


= m^rin-, = \, iuT^ 


0.243588Cbr3 


Pc 


CDM density 


0.243588Ccr3 


Pv 


Neutrino density (massive) 


0.2435885^T3 


7 
Pv 


Neutrino density (massless) 


f (tt)*''%7~ 0.681322p^ 


0.448292r* 


Pm 


Total matter density 


Pb + Pc + P,. = "7? = ^(?b + 5c + 5.)T3 = a;AX 


0.243588Cr3 


Pk 


Curvature density 


3k _ _3_^rp2 


-0.119366kT2 


A 


Expansion factor since equality 


a/acq = Aa^i/s ~ 0.1375145'*/3p;;^/^ 


Si 3.5 X 10^ today 




Dark energy/matter ratio 


PA/Pm = (A/Aa)3 = 


~ 7/3 today 




Hubble reference rate (mere unit) 


H/h = 100km s-^Mpc"! 


(9.7779 Gyr)~i 


Ph 


Hubble reference density (mere unit) SHJ/SttG = 3(100km s-iMpc-i)2/87rG 


1.87882 X 10-26kg/m3 




Baryon density parameter 


^bh^ = Pb/Ph = S.bn^/Ph = (2C(3)/7r2)5bT3 


0.023 ± 0.001 today 


'-'-'cdm 


Cold dark matter density parameter 


^cdrafl^ = Pc/ph = ^cTlj/ph = {2C{3) /n'^ Ph)^,T^ 


0.12 ± 0.01 today 




Neutrino density parameter 


n,h^ = p,/p^ = i^n-,lpn = (2C(3)/7r2ph)C>.T3 


< 0.01 today 


l^d 


Dark matter density parameter 


Odh2 = cc^^dm +i^. = (2C(3)/7r2ph)(§b + ?c)T3 


Si 0.023 ± 0.001 today 




Matter density parameter 


n^h? = c^b + '-'cdm +^u = (2C(3)/7r2p;,)(5b + Cc + S,i.)T^ 0.14 ± 0.01 today 


LOA 


Dark energy density parameter 


Oa/i2 = pA/Ph 


0.34 ± 0.09 




Photon density parameter 


p^/pt^ = (7r2/15ph)T4 Si 1.80618 X 

[^(PA+Pm + P7+Pk)]^^^ 


0.0000247 ± 0.0000004 today 


TJ 

rl 


Hubble parameter 


Si (lOGyr)-i today 


h 


Dimensionless Hubble parameter 




Si 0.7 today 


t 


Age of Universe 


/g"//(a')"Mln a' 


Si 14 Gyr today 


tRD 


Age during radiation era 




Si 0.225492/T2 


tMD 


Age during matter era 


Si 1/[12C(3)]1/251/2t3/2 


Si 0.2633C-i/2y-3/2 


tAD 


Age during vacuum era 


..^ 3 £l/^T 




Qa 


Vacuum density ratio 


i^a/'i^ = PKjh^Ph 


Si 0.7 today 




Matter density ratio 


'^m/h'^ (analogously, Qb = '^b/'^^i f^d = '^d/'^^i eic). 


Si 0.3 today 




Spatial curvature parameter 


1 + K(T//i-)2 


1.01 ± 0.02 today 


A, 


Scalar power normalization 


[(3/27r)(T/lAjK)Q]2/800 


Si 0.8 today 



4. Fecundity: There is no fine-tuning, because intel- 
ligent life of some form will emerge under extremely 
varied circumstances. 

Options 1, 2, and 4 tend to be preferred by physicists, 
with recent developments in inflation and high-energy 
theory giving new popularity to option 2. 

Like relativity theory and quantum mechanics, the 
theory of inflation has not only solved old problems, 
but also widened our intellectual horizons, arguably 
deepening our understanding of the nature of physi- 
cal reality First of all, inflation is generically eternal 
mmiHlliiiSEi, so that even though infla- 
tion has ended in the part of space that we inhabit, it 
still continues elsewhere and will ultimately produce an 
infinite number of post-inflationary volumes as large as 
ours, forming a cosmic fractal of sorts. Second, these re- 
gions may have different physical properties. This can 
occur in fairly conventional contexts involving symmetry 
breaking, without invoking more exotic aspects of eternal 



inflation. (Which, for example, might also allow differ- 
ent values of the axion dark matter density parameter 

I33L Is^l-) More dramatically, a common feature of 
much string theory related model building is that there 
is a "landscape" of solutions, corresponding to spacetime 
conflgurations involving different values of both seem- 
ingly continuous parameters (Table 1) and discrete pa- 
rameters (specifying the spacetime dimensionality, the 
gauge group/particle content, etc.), some or all of which 
may vary across the landscape 35, 36, 37, 38, 39]. 

If correct, eternal inflation might transform that po- 
tentiality into reality, actually creating regions of space 
realizing each of these possibilities. Generically each re- 
gion where inflation has ended is inflnite in size, therefore 
potentially fooling its inhabitants into mistaking initial 
conditions for fundamental laws. Inflation may thus in- 
dicate the same sort of shift in the borderline between 
fundamental and effective laws of physics (at the expense 
of the former) previously seen in theoretical physics, as 
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Effective laws 

("Bylaws", "initial conditions")^ 



Kepler, 
Newton 



* Planets are spherical 

* Orbits are circular 




Initial matter distribution 
(A^i, J^i, ...) 




* SU(3)xSU(2)xU(l) symmetry 

* Dimensionality of space 

* Constants (a, py^, ...) 




TOE equations 



.Fundamental laws. 



FIG. 1: The shifting boundary (horizontal hnes) between fun- 
damental laws and environmental laws/effective laws/initial 
conditions. Whereas Ptolemy and others sought to explain 
roughly spherical planets and circular orbits as fundamental 
laws of nature, Kepler and Newton reclassified such properties 
as initial conditions which we now understand as a combina- 
tion of dynamical mechanisms and selection effects. Classical 
physics removed from the fundamental law category also the 
initial conditions for the electromagnetic field and all other 
forms of matter and energy (responsible for almost all the 
complexity we observe), leaving the fundamental laws quite 
simple. A prospective theory of everything (TOE) incorporat- 
ing a landscape of solutions populated by infiation reclassifies 
important aspects of the remaining "laws" as initial condi- 
tions. Indeed, those laws can differ from one post-inflationary 
region to another, and since inflation generically makes each 
such region enormous, its inhabitants might be fooled into 
misinterpreting regularities holding within their particular re- 
gion as Universal (that is, multiversal) laws. Finally, if the 
Level IV multiverse of all mathematical structures j23 exists, 
then even the "theory of everything" equations that physi- 
cists are seeking are merely local bylaws in Rees' terminology 
[TH . that vary across a wider ensemble. Despite such retreats 
from ab initio explanations of certain phenomena, physics has 
progressed enormously in explanatory power. 



cal physics removed all initial conditions from its predic- 
tive purview, its explanatory power inspired awe. While 
quantum mechanics dashed hopes of predicting when a 
radioactive atom would decay, it provided the founda- 
tions of chemistry, and it predicts a wealth of surprising 
new phenomena, as we continue to discover. 



C. Testing fundamental theories observationally 

Let us group the 31 parameters of Table 1 into a 31- 
dimensional vector p. In a fundamental theory where 
inflation populates a landscape of possibilities, some or 
all of these parameters will vary from place to place as 
described by a 31-dimensional probability distribution 
/(p). Testing this theory observationally corresponds 
to confronting that theoretically predicted distribution 
with the values we observe. Selection effects make this 
challenging |^ if any of the parameters that can 

vary affect the formation of (say) protons, galaxies or ob- 
servers, then the parameter probability distribution dif- 
fers depending on whether it is computed at a random 
point, a random proton, a random galaxy or a random 
observer 0, Q| . A standard application of conditional 
probabilities predicts the observed distribution 



/(P) « /prior (p)/scloc(p), 



(1) 



where /prior (p) is the theoretically predicted distribution 
at a random point at the end of inflation and /seiec(p) is 
the probability of our observation being made at that 
point. This second factor /scioc(p), incorporating the 
selection effect, is simply proportional to the expected 
number density of reference objects formed (say, protons, 
galaxies or observers). 

Including selection effects when comparing theory 
against observation is no more optional than the correct 
use of logic. Ignoring the second term in equation can 
even reverse the verdict as to whether a theory is ruled 
out or consistent with observation. On the other hand, 
anthropic arguments that ignore the first term in equa- 
tion are likewise spurious; it is crucial to know which 
of the parameters can vary, how these variations are cor- 
related, and whether the typical variations are larger or 
smaller than constraints arising from the selection effects 
in the second term. 



illustrated in Figure ^ 

There is quite a lot that physicists might once have 
hoped to derive from fundamental principles, for which 
that hope now seems naive and misguided 24]. Yet it 
is important to bear in mind that these philosophical 
retreats have gone hand in hand with massive progress 
in predictive power. While Kepler and Newton discred- 
ited ah initio attempts to explain planetary orbits and 
shapes with circles and spheres being "perfect shapes", 
Kepler enabled precise predictions of planetary positions, 
and Newton provided a dynamical explanation of the ap- 
proximate sphericity of planets and stars. While classi- 



D. A case study: cosmology and dark matter 

Examples where we can compute both terms in equa- 
tion are hard to come by. Predictions of fundamen- 
tal theory for the first term, insofar as they are plau- 
sibly formulated at present, tend to take the form of 
functional constraints among the parameters. Famil- 
iar examples are the constraints among couplings aris- 
ing from gauge symmetry unification and the constraint 
^QCD ~ arising from Peccei-Quinn symmetry. At- 
tempts to predict the distribution of inflation-related 
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cosmological parameters are marred at present by reg- 
ularization issues related to comparin g in finite volumes 

PEEllillillillilliilillilliiM Addi- 
tional difficulties arise from our limited understanding as 
to what to count as an observer, when we consider varia- 
tion in parameters that affect the evolution of life, such as 
(jTip, a, P), which approximately determine all properties 
of chemistry. 

In this paper, we will focus on a rare example where 
where there are no problems of principle in computing 
both terms: that of cosmology and dark matter, involv- 
ing variation in the parameters (^c, pA, Q) from Table 1, 
i.e., the dark matter density parameter, the dark energy 
density and the seed fluctuation amplitude. Since none of 
these three parameters affect the evolution of life at the 
level of biochemistry, the only selection effects we need to 
consider are astrophysical ones related to the formation 
of dark matter halos, galaxies and stable solar systems. 
Moreover, as discussed in the next section, we have spe- 
cific well-motivated prior distributions for pA and (for 
the case of axion dark matter) ^c- Making detailed dark 
matter predictions is interesting and timely given the ma- 

t efforts underway to detect dark matter both directly 
and indirectly [53| and the prospects of discovering 
supersymmetry and a WIMP dark matter candidate in 
Large Hadron Collider operations from 2008. 

For simplicity, we do not include any of the cur- 
rently optional cosmological parameters, i.e., we take 
ris = fitot = 1, a„ = r = nt = 0, w = -1. The 
remaining two non-optional cosmological parameters in 
Table 1 are the density parameters for neutrinos and 
for baryons. It would be fairly straightforward to 
generalize our treatment below to include along the 
lines of [s^, ; since it too affects /seiec only through 
astrophysics and not through subtleties related to bio- 
chemistry. Here, for simplicity, we will ignore it; in any 
case, it has been observed to be rather unimportant cos- 
mologically {^u <^ ^c)- When computing cosmological 
fluctuation growth, we will also make the simplifying ap- 
proximation that ^ Cc (so that ^ ~ ^c), although 
we will include as a free parameter when discussing 
galaxy formation and solar system constraints. (For very 
larg e £b , structure formation can change qualitatively; 
see [103 •) We will see below that ic/ih 3> 1 is not only 
observationally indicated (Table 1 gives fc/^b ~ 6), but 
also emerges as the theoretically most interesting regime 
if ^b is considered fixed. 

The rest of this paper is organized as follows. In Sec- 
tion ^ we discuss theoretical predictions for the first 
term of equation the prior distribution /prior- In Sec- 
tion lllll we discuss the second term /sciec, computing the 
selection effects corresponding to halo formation, galaxy 
formation and solar system stability. We combine these 
results and make predictions for dark-matter-related pa- 
rameters in Section IIVI summarizing our conclusions in 
Section A number of technical details are relegated 
to Appendix A. 



II. PRIORS 

In this section, we will discuss the first term in equa- 
tion QJ, specifically how the function /prior (p) depends 
on the parameters fc, PA and Q. In the case of fc, we will 
consider two dark matter candidates, axions and WIMPs. 



A. Axions 

The axion dark matter model offers an elegant example 
where the prior probability distribution of a parameter 
(in this case ^c) can be computed analytically. 

The strong CP problem is the fact that the dimension- 
less parameter 0qcd in Table 1, which parameterizes a 
potential CP-violating term in quantum chromodynam- 
ics (QCD), the theory of the strong interaction, is so 
small. Within the standard model, O^cd is a periodic 
variable whose possible values run from to 27r, so its 
natural scale is of order unity. Selection effects are of 
little help here, since values of 6'qcd far larger than the 
observed bound |6'qcd| ^ 10~^ would seem to have no 
serious impact on life. 

Peccei and Quinn [ssj ] introduced microphysical models 
that address the strong CP problem. Their models ex- 
tend the standard model so as to support an appropriate 
(anomalous, spontaneously broken) symmetry. The sym- 
metry is called Peccei-Quinn (PQ) symmetry, and the en- 
ergy scale at which it breaks is called the Peccei-Quinn 
scale. Weinberg 5Gj and Wilczek independently re- 
alized that Peccei-Quinn symmetry implies the existence 
of a field whose quanta are extremely light, extremely 
feebly interacting particles, known as axions. Later it 
was shown that axions provide an interesting dark mat- 
ter candidate 0,0,0. 

Major aspects of axion physics can be understood by 
reference to a truncated toy model where O^cd is the com- 
plex phase angle of a complex scalar field $ that develops 
a potential of the type 

y($) = - flY + k%^{l - /-iRe<I>), (2) 

where Aqcd ~ 200 MeV is ultimately determined by the 
parameters in Table 1; roughly speaking, it is the energy 
scale where the strong coupling constant as (Aqcd) = 1- 

At the Peccei-Quinn (PQ) symmetry breaking scale 
fa, assumed to be much larger than Aqcd, this complex 
scalar field $ feels a Mexican hat potential and seeks to 
settle toward a minimum ($) = /ae'^i'^'*. In the context 
of cosmology, this will occur at temperatures not much 
below fa. Initially the angle Oqcd = 6*0 is of negligible 
energetic significance, and so it is effectively a random 
field on superhorizon scales. The angular part of this 
field is called the axion field a — fadqcd- As the cosmic 
expansion cools our universe to much lower temperatures 
approaching the QCD scale, the approximate azimuthal 
symmetry of the Mexican hat is broken by the emergence 
of the second term, a periodic potential 1 — /~^Re$ = 
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1 — cos Oqcd = 2 sin^ 9qcd (induced by QCD instantons) 
whose minimum corresponds to no strong CP-violation, 
i.e., to 6'qcd = 0. The axion field oscillates around this 
minimum like an underdamped harmonic oscillator with 
a frequency nia corresponding to the second derivative of 
the potential at the minimum, gradually settling towards 
this minimum as the oscillation amplitude is damped by 
Hubble friction. That oscillating field can be interpreted 
as a Bose condensate of axions. It obeys the equation of 
state of a low-pressure gas, which is to say it provides a 
form of cold dark matter. 

By today, 6'qcd is expected to have settled to an an- 
gle within about 10~^® of its minimum [63, comfortably 
below the observational limit \6qcd\ ^ 10~^, and thus dy- 
namically solving the strong CP problem. (The exact 
location of the minimum is model-dependent, and not 
quite at zero, but comfortably small in realistic models 

The axion dark matter density per photon in the cur- 
rent epoch is estimated to be |58l l5£ 



i e • 2 Co 



C* ~ /a ■ 



(3) 



If axions constitute the cold dark matter and the Peccei- 
Quinn phase transition occurred well before the end of 
inflation, then the measurement ~ 3 x 10~^* thus im- 
plies that 

-1/2 / . ^ -1/2 



fa - 



sm ■ 



10'^ GeV X sin 



(4) 

where is the the initial misalignment angle of the axion 
field in our particular Hubble volume. 

Frequently it has been argued that this implies fa ^ 
10^^ GeV, ruHng out GUT scale axions with fa ^ 
10^^ GeV. Indeed, in a conventional cosmology the hori- 
zon size at the Peccei-Quinn transition corresponds to 
a small volume of the universe today, and the observed 
universe on cosmological scales would fully sample the 
random distribution 6q. However, the alternative possi- 
bility that 1 00 1 <C 1 over our entire observable universe 
was pointed out already in 58] . It can occur if an epoch 
of inflation intervened between Peccei-Quinn symmetry 
breaking and the present; in that case the observed uni- 
verse arises from within a single horizon volume at the 
Peccei-Quinn scale, and thus plausibly lies within a cor- 
relation volume. Linde |32j| argued that if there were 
an anthropic selection effect against very dense galaxies, 
then models with fa > 10^^ GeV and |6io| < 1 might 
indeed be perfectly reasonable. Several additional as- 
pects of this scenario were discussed in |33i, ,3^ . Much 
of the remainder of this paper arose as an attempt to 
better ground its astrophysical foundations, but most of 
our considerations are of much broader application. 

We now compute the axion prior /prior (Cc)- Since the 
symmetry breaking is uncorrelated between causally dis- 
connected regions, Oq is for all practical purposes a ran- 
dom variable that varies with a uniform distribution be- 
tween widely separated Hubble volumes. Without loss of 



generality, we can take the interval over which varies 
to be < 00 ^ I"- This means that the probability of ^ 
being lower than some given value ^o < is 

< Co) = Pi^* sin^ ^ < = P (^^0 < 2 sin-i 



2 ^^/^ 



(5) 



Differentiating this expression with respect to ^o gives the 
prior probability distribution for the dark matter density 



fdO 



1 



(6) 



For the case at hand, we only care about the tail of the 
prior corresponding to unusually small 9o, i.e., the case 
S, <^ for which the probability distribution reduces to 
simply 



fdO « 



1 



(7) 



Although this may appear to favor low ^, the probability 
per logarithmic interval cx ^-nd it is obvious from 
equation (0 that the bulk of the probability lies near the 
very high value ^ ~ ^* . 

A striking and useful property of equation ||7J) is that it 
contains no free parameters whatsoever. In other words, 
this axion dark matter model makes an unambiguous pre- 
diction for the prior distribution of one of our 31 pa- 
rameters, ^c- Since the axion density is negligible at 
the time of inflation, this prior is immune to the infla- 
tionary measure-related problems discussed in 0, and 
no inflation-related effects should correlate with other 
observable parameters. Moreover, this conclusion applies 
for quite general axion scenarios, not merely for our toy 
model — the only property of the potential used to derive 
the ^^^-scaling is its parabolic shape near any mini- 
mum. Although many theoretical subtleties arise regard- 
ing the axion dark matter scenario in the contexts of in- 
flation, supersymmetry and string theory 62, 63, 64, 65|, 

the prior appears as a robust consequence of the 

hypothesis fa > lO^^ Q^y 

We conclude this section with a brief discussion of 
bounds from axion fluctuations. Like any other massless 
field, the axion field a = faO acquires fluctuations of or- 
der H during inflation, where H ^ -^inf/"^pi ~ ^fnf 
Eini is the inflationary energy scale, so SOq ~ ^ini/fa- 
For our j^ol ^ 1 case, equation (O gives Bq ~ Cc^^/Za, 
so we obtain the axion density fluctuations amplitude 



sc 



2500 



(8) 



Such axion isocurvature fluctuations (see, e.g., for 
a review) would contribute acoustic peaks in the cos- 
mic microwave background (CMB) out of phase with 
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the those from standard adiabatic fluctuations, allow- 
ing an observational upper bound Qa ^ O.SQ ~ 10"^ 
to be placed E^- Combining this with the ob- 
served ^c-value from Table 1 gives the bound Ef^^fa ^ 
Qa^y^ ^ 10~^^, bounding the inflation scale. The tra- 
ditional value fa ^ 10^^ GeV ^ 10 ~J g ives the famil- 
iar bound Einf < 10-s - 10" GeV A higher 
fa gives a tighter limit on the inflation scale: increas- 
ing fa to the Planck scale {fa ^ 1) lowers the bound to 
^'inf ^ 10~^ ~ 10^*^ GeV — the constraint grows stronger 
because the denominator in 59q/9q must be smaller to 
avoid an excessive axion density. 

For comparison, inflationary gravitational waves have 
amplitude rQ ^ H ^ -^inf' ^'^ they are unobservably 
small unless E'inf ^ 10^^ GeV. Although various loop- 
holes to the axion fluctuation bounds have been proposed 
(see e.g., [HE E3, llS it is interesting to note that 
the simplest axion dark matter models therefore make 
the falsifiable prediction that future CMB experiments 
will detect no gravitational wave signal • 



B. WIMPs 

Another popular dark matter candidate is a stable 
weakly interacting massive particle (WIMP), thermally 
produced in the early universe and with its relic abun- 
dance set by a standard freezeout calculation. See, e.g., 
[Z^l73| for reviews. Stability could, for instance, be en- 
sured by the WIMP being the lightest supersymmetric 
particle. 

At relevant (not too high) temperatures the thermally 
averaged WIMP annihilation cross section takes the form 
III 



7- 



(9) 



where 7 is a dimensionless constant of order unity. In 
our notation, the WIMP number density is Jiwimp = 
Pwimp/TOwimp = "7Cwimp/m-wimp- The WIMP freczout is 
determined by equating the WIMP annihilation rate P = 
(crw) riwimp = (o'u)n^^wimp/"iwimp with the radiation- 
dominated Hubble expansion rate H = (87rp^/3)^/^. 
Solving this equation for ^wimp and substituting the ex- 
pressions for and p~f from Tables 2 and 3 gives 



297ri 
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(10) 



The WIMP freezeout temperature is typically found to 
be of order T ~ TOwimp/20 [l^l- If we further assume 
that the WIMP mass is of order the electroweak scale 
(wwimp v) and that the annihilation cross section pref- 
actor 7 ~ 1, then equation H10() gives 



for the measured values of v and g from Table 1. This 
well-known fact that the predicted WIMP abundance 
agrees qualitatively with the measured dark matter den- 
sity is a key reason for the popularity of WIMP dark 
matter. 

In contrast to the above-mentioned axion scenario, we 
have no compelling prior for the WIMP dark matter den- 
sity parameter ^wimp- Let us, however, briefly explore the 
interesting scenario advocated by, e.g., |7J|, where the 
theory prior determines all relevant standard model pa- 
rameters except the Higgs vacuum expectation value v, 
which has a broad prior distribution.'^ (It should be said, 
however, that the connection Wwimp ~ f is somewhat ar- 
tiflcial in this context.) It has recently been shown that v 
is subject to quite strong microphysical selection effects 
that have nothing to do with dark matter, as nicely re- 
viewed in . As pointed out by IZ3 j changing v 
up or down from its observed value vq by a large factor 
would correspond to a dramatically less complex universe 
because the slight neutron-proton mass difference has a 
quark mass contribution (m^ — m„) cx v that slightly 
exceeds the extra Coulomb repulsion contribution to the 
proton mass: 

1. For v/va < 0.5, protons (uud) decay into neutrons 
(udd) giving a universe with no atoms. 

2. For v/vq > 5, neutrons decay into protons even in- 
side nuclei, giving a universe with no atoms except 
hydrogen. 

3. For v/vq > 10'^, protons decay into A"''^ (uuu) par- 
ticles, giving a universe with only Helium-like 
atoms. 

Even smaller shifts would qualitatively alter the synthe- 
sis of heavy elements: For v/vq < 0.8, diprotons and 
dineutrons are bound, producing a universe devoid of, 
e.g., hydrogen. For v/vq > 2, deuterium is unstable, 
drastically altering standard stellar nucleosynthesis. 

Much stronger selection effects appear to result from 
carbon and oxygen production in stars. Revisiting the 
issue first identified by Hoyle [78| with numerical nu- 
clear physics and stellar nucleosynthesis calculations, [t^ 
quantified how changing the strength of the nucleon- 
nucleon interaction altered the yield of carbon and oxy- 
gen in various types of stars. Combining their results 
with those of [801 that relate the relevant nuclear physics 
parameters to v gives the following striking results: 

1. For v/vq < 0.99, orders of magnitude less carbon is 
produced. 



10^ 



lO'^ 



10 



-28 



(11) 



^ If only two microphysical parameters from Table 1 vary by many 
orders of magnitude across an ensemble and are anthropically 
selected, one might be tempted to guess that they are pA ~ 
since they differ most dramatically 
from unity. A broad prior for — /x^ would translate into a broad 
prior for v that could potentially provide an anthropic solution 
to the so-called "hierarchy problem" that v ~ 10"^^ <^ 1. 



2. For v/vq > 1.01, orders of magnitude less oxygen 
is produced. 

Combining this with equation I|1HI , we see that this could 
potentially translate into a percent level selection effect 
on ^wimp oc v'^. 

In the above-mentioned scenario where v has a broad 
prior whereas the other particle physics parameters (in 
particular g) do not [t^, the fact that microphysical se- 
lection effects on v are so sharp translates into a narrow 
probability distribution for ^wimp via equation (|ll|l . As 
we will see, the astrophysical selection effects on the dark 
matter density parameter are much less stringent. We 
should emphasize again, however, that this constraint 
relies on the assumption of a tight connection between 
Cwimp and w, which could be called into question. 



C. PA and Q 

As d iscussed in de t ail i n the hterature {e.g., 0, 0, 
H m Eol Eol Eoi [ll3), there are plausible reasons 
to adopt a prior on pA that is essentially constant and 
independent of other parameters across the narrow range 
where |pa| < Q^S,'*^ ~ 10^^^^ where /scicc is non-neghgible 
(see Section HlB . The conventional wisdom is that since 
Pa is the difference between two much larger quantities, 
and Pa = has no evident microphysical significance, no 
ultrasharp features appear in the probability distribution 
for PA within 10"^^'^ of zero. That argument holds even 
if Pa varies discretely rather than continuously, so long 
as it takes ^ 10^^'^ different values across the ensemble. 

In contrast, calculations of the prior distribution for Q 
from inflation are fraught with considerable uncertainty 
[6, 81, 82] . We therefore avoid making assumptions about 
this function in our calculations. 
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FIG. 2: Many selection effects that we discuss arc conveniently 
summarized in the plane tracking the virial temperatures and den- 
sities in dark matter halos. The gas cooling requirement prevents 
halos below the heavy black curve from forming galaxies. Close en- 
counters make stable solar systems unlikely above the downward- 
sloping line. Other dangers include collapse into black holes and 
disruption of galaxies by supernova explosions. The location and 
shape of the small remaining region (unshaded) is independent of 
all cosmological parameters except the baryon fraction. 



III. SELECTION EFFECTS 

We now consider selection effects, by choosing our "se- 
lection object" to be a stable solar system, and focusing 
on requirements for creating these. In line with the pre- 
ceding discussion, our main interest will be to explore 
constraints in the 4-dimensional cosmological parameter 
space (ifb, i^c, PA, Q). Since we can only plot one or two 
dimensions at a time, our discussion will be summarized 
by a table (Table 4) and a series of figures showing var- 
ious 1- and 2-dimensional projections: (Cb,Cc), (PA,QO: 

Many of the the physical effects that lead to these con- 
straints are summarized in Figure |21 and Figure |3| show- 
ing temperatures and densities of galactic halos. The 
constraints in this plane from galaxy formation and solar 
system stability depend only on the microphysical pa- 
rameters (mp,Q;,/?) and sometimes on the baryon frac- 
tion ^b/Ccj whereas the banana-shaped constraints from 
dark matter halo formation depend on the cosmologi- 
cal parameters (^b, Cc, £,vt PKtQ), so combining them con- 
strains certain parameter combinations. Crudely speak- 
ing, /soioc(p) will be non- negligible only if the cosmolog- 
ical parameters are such that part of the "banana" falls 
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FIG. 3: Same as the previous figure, but including the banana- 
shaped contours showing the the halo formation/destruction rate. 
Solid/greenish contours correspond to positive net rates (halo 
production) at 0.5, 0.3, 0,2, 0.1, 0.03 and 0.01 of maximum, 
whereas dashed/reddish contours correspond to negative net rates 
(halo destruction from merging) at -0.3, -0,2, -0.1, -0.03 and - 
0.01 of maximum, respectively. The heavy black contour corre- 
sponds to zero net formation rate. The cosmological parameters 
(?bi ?c, ^i/, , Q, pa) have a strong effect on this banana: gi, and 
shift this banana-shape vertically, Q shifts if along the parallel di- 
agonal lines and p\ cuts it off below 16pA (see the next figure). 
Roughly speaking, there are stable habitable planetary systems 
only for cosmological parameters where a greenish part of the ba- 
nana falls within the allowed white region from Figure 1^ 

within the "observer-friendly" (unshaded) region sand- 
wiched between the galaxy formation and solar system 
stability constraints 

In the following three subsections, we will now dis- 
cuss the three above-mentioned levels of structure for- 
mation in turn: halo formation, galaxy formation and 
solar-system stability. 



A. Halo formation and the distribution of halo 
properties 

Previous studies {e.g., |1 E H HI [lol [loi HM 

have computed the total mass fraction collapsed 
into halos, as a function of cosmological parameters. 
Here, however, we wish to apply selection effects based on 
halo properties such as density and temperature, and will 
compute the formation rate of halos (the banana-shaped 
function illustrated in Figure|31) in terms of dimensionless 
parameters alone. 
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FIG. 4: Same as Figure |21 but showing effect of changing the 
cosmological parameters. Increasing the matter density parameter 
5 shifts the banana up by a factor (top left panel), increasing 
the fluctuation level Q shifts the banana up by a factor and 
to the right by a factor Q (top right panel), increasing the baryon 
fraction = 5b /C shifts the banana up by a factor and affects 
the cooling and 2nd generation constraints (bottom left panel), and 
increasing the dark energy density pA bites off the banana below 
Pvir = 16pA (bottom right panel). In these these four panels, the 
parameters have been scaled relative to their observed values as 
follows: up by 1/4, 2/4,...,9/4 orders of magnitude (5), up by 1/3, 
2/3,...,9/3 orders of magnitude (Q), down by 1 and 2 orders of 
magnitude (/b), and down by 0, 7, 12 and oo orders of magnitude 
(pa; other panels have pA = 0). 



1. The dependence on mass and time 

As previously explained, we will focus on the dark- 
matter-dominated case ^ Cbi ^ Ci^, so we ignore 
massive neutrinos and have ^ = + Cc ~ ■ For our cal- 
culations, it is convenient to define a new dimensionless 
time variable 



_ PA 



1 



and a new dimensionless mass variable 



(12) 



(13) 



In terms of the usual cosmological scale factor a, our new 
time variable therefore scales as a; oc a"^. It equals unity 
at the vacuum domination epoch when linear fluctuation 
growth grinds to a halt. The horizon mass at matter- 
radiation equality is of order [s^, so fi can be in- 
terpreted as the mass relative to this scale. It is a key 
physical scale in our problem. It marks the well-known 



break in the matter power spectrum; fluctuation modes 
on smaller scales entered the horizon during radiation 
domination, when they could not grow. 

We estimate the fraction of matter collapsed into dark 
matter halos of mass M > by time x using the 

standard Press-Schechter formalism |8q . which gives 



F{fi,x) = erfc 



Scix) 



(14) 



Here (t(/z, x) is the r.m.s. fluctuation amplitude at time 
a; in a sphere containing mass M = so F is the 

probability that a fluctuation lies 5c standard deviations 
out in the tail of a Gaussian distribution. As shown in 
Appendix lA 21 a is well approximated as^ 



(j{n,x) 



' 1/3 

Pa 



where 



45 • 21/3^(3)4/3 



ri4/3 



1+f (^) 



4/3 



0.206271 



(15) 



(16) 



and the known dimensionless functions s(/x) and G/^{x) 
do not depend on any physical parameters, s(^) and 
G\{x) give the dependence on scale and time, respec- 
tively, and appear in equations (|A13(I and ljAl|l in Ap- 
pendix A. The scale dependence is s(/i) ~ fi^^^^ on large 
scales n ^ 1, saturating to only logarithmic growth 
towards small scales for /i <C 1. Fluctuations grow 



as Ga(x) 



,1/3 



cx a for X <^ 1 and then asymp- 



tote to a constant amplitude corresponding to Goo = 
5r (I) r (I) /3V^ w 1.43728 as a ^ oo and dark energy 
dominates. (This is all for pA > 0; we will treat pA < 
in Section IIII A 31 and find that our results are roughly 
independent of the sign of p\, so that we can sensibly 
replace pA by \pa\-) 

Returning to equation the collapse density 

threshold Sc{x) is defined as the linear perturbation 
theory overdensity that a top-hat-averaged fluctuation 
would have had at the time x when it collapses. It was 
computed numerically in j^, and found to vary only 
very weakly (by about 3%) with time, dropping from the 
famihar cold dark matter value Sc{0) = (3/20)(127r)2/3 « 
1.68647 early on to the limit ^^(oo) = {9/5)2-^^^G^ « 
1.62978 in the infinite future. Here we simply ap- 
proximate it by the latter value: 



Sc{x) « Sc{oo) 



9T 



3^1/222/3 



1.62978. 



(17) 



4 Although the baryon density affects a mainly via the sum ^ = 
Cb + Cc, there is a slight correction because fluctuations in the 
baryon component do not grow between matter-radia tion equal- 
ity and the drag epoch shortly after recombination |87j . Since 
the resulting correction to the fluctuation growth factor (~ 15% 
for the observed baryon fraction 5b /5c ~ 1/6) is negligible for 
our purposes, we ignore it here. 
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Substituting equations H15|) and H17() into equation (|14() . 
we thus obtain the coUapsed fraction 



F(/Lt, x) = erfc 



(18) 



0.8 



where 



^ = = -^^ , . .^r — ^ w 5.58694. 



30\/2C(3)4/3 



(19) I 



Below we will occasionally find it useful to rewrite equa- ^ 
tion l(TH|l as 



F(^, x) = erfc 



Ga{x)s{p) 



4n3 



(20) 



Let us build some intuition for equation p8|l . It tells us 
that early on when x ~ 0, no halos have formed (F « 0), 
and that as time passes and x increases, small halos form 
before any large ones since s(/i) is a decreasing function. 
Moreover, we see that since Ga{x) and s{p) are at most 
of order unity, no halos will ever form if pA ^ £,'^Q^- We 
recognize the combination p^ = ^^Q^ as the character- 
istic density of the universe when halos would form in 
the absence of dark energy js^. If pa S> p*, then dark 
energy dominates long before this epoch and fluctuations 
never go nonlinear. Figure illustrates that the density 
distributions corresponding to equation (|2U|I are broadly 
peaked around pvir ^ lO^p, when ^ pA, are exponen- 
tially suppressed for p* <ti pa, and in all cases give no 
halos with pvir < 16pA. 



2. The dependence on temperature and density 

We now discuss how our halo mass p and formation 
time X transform into the astrophysically relevant pa- 
rameters (Tvir,n.vir) appearing in Figure|21 

As shown in Appendix lA 31 halos that virialize at time 
X have a characteristic density of order 
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FIG. 5: Halo density distribution for various values of pt/pi\, 
where p* = [s(/x)/A]^p, ~ Q^?'*. The dashed curves show the cu- 
mulative distribution F = erfc [(pA//3*)"'^''^/GA(a;[pvir])] and the 
solid curves show the probability distribution — 9F/91gpvir for 
(peaking from right to left) p*/pA = 10^, 10"', 10^, 10^, 10, 1 
(heavy curves), 0.1 and 0.05. Note that no halos with pvir < 16pA 
are formed. 



Inverting equations H21|) and (|22|l gives 



/ £4^3. 
S vir 

"l^Pvir 



97r^ 



( Pv 



\IQPA 



(23) 
(24) 



For our applications, the initial gas temperature will be 
negligible and the gas density will trace the dark mat- 
ter density, so until cooling becomes important (Sec- 
tion the proton number density is simply 



~ Pvir 



(25) 



i.e., essentially the larger of the two terms 16pA and 
187r^Pi„(a;). For a halo of total (baryonic and dark 
matter) mass M in Planck units, this corresponds to 
a characteristic size R ~ (M/pvir)^^'^, velocity i^vir ~ 
{M/Ry/^ ~ (M^pvir)^/^ and virial temperature Tvir ~ 



T ■ 

^ vir 



/ 16pAM^ 



so 



1/3 



9^ 
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(22) 



According to the Press-Schechter approximation, the 
derivative — g^^^ (p, x) can be interpreted as the so-called 
mass function, i.e., as the distribution of halo masses at 
time X. We can therefore interpret the second derivative 



Ux{P,x) 



d{\gp)d{\gx) 



(26) 



as the net formation rate of halos as a function of mass 
and time. Transforming this rate from (lgp,lgx)-space 
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to (lgTvir,lgPvii)-space, we obtain the function whose 
banana-shaped contours are shown in Figure 13 



\J\Ux{i^,x) 



\J\ 



a(ig^)a(igx)' 



(27) 



where J is a Jacobian determinant that can be ignored 
safely.^ 

Let us now build some intuition for this important 
function f {T^ir , n^ii) ■ First of all, since equation lpT|l 
shows that pvir c>c nvir decreases with time and is jj,- 
independent, we can reinterpret the vertical rivir-axis in 
Figure 01 as simply the time axis: as our Universe ex- 
pands, halos can form at lower densities further and fur- 
ther down in the plot. Since nothing ever forms with 
Pvir < 16/9A, the horizontal line rivir = ISpA^b/C^^p cor- 
responds to i — > oo. Second, / is the net formation rate, 
which means that it is negative if the rate of formation 
of new halos of this mass is smaller than the rate of de- 
struction from merging into larger halos. The destruction 
stems from the fact that, to avoid double-counting, the 
Press-Schechter approximation counts a given proton as 
belonging to at most one halo at a given time, defined 
as the largest nonlinear structure that it is part of. ^ 
Figure IHl shows that halos of any given mass (defined by 
the lines of slope 3) are typically destroyed in this fash- 
ion some time after its formation unless pA-domination 
terminates the process of fluctuation growth. 

We will work out the detailed dependence of this ba- 
nana on physical parameters in the next section. For 
now, we merely note that Figure |31 shows that the first 
halos form with characteristic density pvir ~ p* ~ S.'^Q^i 
n-vii ^ ^h£,^Q^ /mp, and (assuming pA ^ p*) the largest 
halos approach virial velocities Wvir ~ Q^^^c and temper- 
atures Tvir ~ Q X TOpC^ [85l |. 



It makes sense to treat ffj.x as a distribution, since 
/ / ffj,xd{ln ii) d(lnx) equals the total collapsed fraction. When 
transforming it, we therefore factor in the Jacobian of the trans- 
formation from (lgfl,lgx) to (Ig Tvir, Ig Pvir)i 



/ 8 Ig M 8 Ig M 

J = aigT eigp 

I dlgx dig X 

\ aigT sigp' 

with determinant 



J=|J| 
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(28) 



(29) 



So the Jacobian is an irrelevant constant J = —3/2 for Pvir 3> 
16pA- Since the entire function / vanishes for pvir < 16pA, 
the Jacobian only matters near the Pvir = 16pA boundary, 
where it has a rather unimportant effect (the divergence is inte- 
grable). Equation I28i shows that away from that boundary, the 
(Ig Tvir, Ig Pvir )-t)anana is simply a linear transformation of the 
(Ig /i, Ig a;)-banana F/dlg fidlgx, with slanting parallel lines of 
slope 3 in Figure |21 corresponding to constant p- values. 
Our treatment could be improved by modeling halo substructure 
survival, since subhalos that harbor stable solar systems may be 
counted as part of an (apparently inhospitable) larger halo. 



3. How TUp, ^b, ?c, Q and pA affect "tfie banana" 

Substituting the preceding equations into equa- 
tion (|27|) gives the explicit expression 

/(Tvir, "^vir; rUp, ^b, Cc, Q, Pa) = 
3a\GA(xfs{pf-2a^]xGA'{x)ns'{tJ') 



v/^FGA(a;)'is(M)*cxp ( 



( _P^ArX_\ 



(30) 



where fi and x are determined by Tvir and nvir via equa- 
tions H23|) , H24|) and (|25|) . This is not very illuminating. 
We will now see how this complicated-looking function / 
of seven variables can be well approximated and under- 
stood as a fixed banana-shaped function of merely two 
variables, which gets translated around by variation of 
"T'pj ^b, Cc and Q, and truncated from below at a value 
determined by pA- 

Let us first consider the limit pA ^ where dark en- 
ergy is negligible. Then a; ^ 1 so that Ga{x) = x^^"^ = 
{pk/ PmY^^ , causing equation ((TH|l to simphfy to 



T(p, x) = erfc 



Ap 



1/3 



^4/3Q,(^) 



(31) 



In the same limit, equation (|21|l reduces to pvir = 
ISTT^Pin, so using equation gives pm = pvir/187r^ = 
^TOpTivir/lSTr^^b- Using this and equation H23|l . we can 
rewrite equation H31|l in the form 



F — erfc 



A 



1/3 



= B 



(18^2)1/3^ 
TTZpTT-vir ^vir 



^1/2- 



Cb^Q^ ' mpQ 



(32) 



where the "standard banana" is characterized by a func- 
tion of only two variables. 



B{X, Y) = erfc 



AX^/^ 



(187r2)i/S [y3/2x-i/2] 



(33) 



Thus B{X, Y) determines the banana shape, and the pa- 
rameters TOp, ^b, £, and Q merely shift it on a log- log 
plot: increasing nip shifts it down and to the right along 
lines of slope —1, increasing ^b shifts it upward, increas- 
ing ^ shifts it upward three times faster, and increasing 
Q shifts it upward and to the right along lines of slope 3. 
The halo formation rate defined by its derivatives (equa- 
tion and plotted in Figure 13 clearly scales in exactly 
the same way with these four parameters, since in the 
Pa — ^ limit that we are considering, the Jacobian J is 
simply a constant matrix. 

Let us now turn to the general case of arbitrary pA. As 
as discussed above, time runs downward in Figure|21since 
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the cosmic expansion gradually dilutes the matter den- 
sity Pin- The matter density completely dwarfs the dark 
energy density at very early times. They key point is that 
since pA has no effect until our Universe has expanded 
enough for the matter density to drop near the dark en- 
ergy density, the part of the banana in Figure El that 
lies well above the vacuum density pA will be completely 
independent of pA- As the matter density pm drops be- 
low Pa {i.e., as x grows past unity), fluctuation growth 
gradually stops. This translates into a firm cutoff below 
Pvir = 16pA in Figure 13 (i.e., below rivir = 16^bPA/"^pC)j 
since equation (|21|) shows that no halos ever form with 
densities below that value. 

Viewed at sensible resolution on our logarithmic plot, 
spanning many orders of magnitude in density, the transi- 
tion from weakly perturbing the banana to biting it off is 
quite abrupt, occurring as the density changes by a factor 
of a few. For the purposes of this paper, it is appropriate 
to approximate the effect of pA as simply truncating the 
banana below the cutoff density. 

Putting it together, we can approximate the non- 
intuitive equation H30() by the much simpler 



/(Tvii., rivir; mp, ^b, Cc, Q, pa) 



16CbPA 

nip 5 



(34) 



where 6 is the Heaviside step function {0{x) = for x < 0, 
6{x) = 1 for a: > 0) and b is the differentiated "standard 
banana" function 



biX,Y) 



92 



igxigy 



B{X,Y). 



(35) 



This is useful for understanding constraints on the cos- 
mological parameters: the seemingly complicated depen- 
dence of the halo distribution on seven variables from 
equation H3()|l can be intuitively understood as the stan- 
dard banana shape from Figure|31being rigidly translated 
by the four parameters nip, ^b, S, and Q and truncated 
from below with a cutoff n^ir > 16^bPA/'Tip^. All this is 
illustrated in Figure 21 which also confirms numerically 
that the approximation of equation ()35|l is quite accurate. 



4- The case pA < 

Above we assumed that pA > 0, but for the purposes 
of this paper, we can obtain a useful approximate gen- 
eralization of the results by simply replacing pA by |pa| 
in equation (|34|l . This is because pA has a negligible ef- 
fect early on when |pa| ^ Pm and a strongly detrimental 
selection effect when |pa| ~ Pm, either by suppressing 
galaxy formation (for pA > 0) or by recoUapsing our uni- 
verse (for PA < 0), in either case causing a rather sharp 
lower cutoff of the banana. As pointed out by Wein- 
berg [s^, one preferentially expects pA > as observed 
because the constraints tend to be slightly stronger for 
negative pA. If pa > 0, observers have time to evolve 



long after |pa| = Pm as long as galaxies had time to 
form before while pA was still subdominant. If pa < 0, 
however, both galaxy formation and observer evolution 
must be completed before pA dominates and recoUapses 
the universe; thus, for example, increasing Q so as to 
make structure form earlier will not significantly improve 
prospects for observers. 



B. Galaxy formation 

1. Cooling and disk formation 

Above we derived the time-dependent (or equivalently, 
density-dependent) fraction of matter collapsed into ha- 
los above a given mass (or temperature), derived from 
this the formation rate of halos of a given temperature 
and density, and discussed how both functions depended 
on cosmological parameters. For the gas in such a halo to 
be able to contract and form a galaxy, it must b e abl e to 
dissipate energy by coohng 96, 97, 98, 99] — see |10(]| for 
a recent review. In the shaded region marked "No cool- 
ing" in Figure 13 the cooling timescale T/T exceeds the 
Hubble timescale H~^J We have computed this familiar 
cooling curve as in [s^ with updated molecular cooling 
from Tom Abel's code based on llOll. which is available 



at Ihttp: //www. tomabel.com/PGas/, From left to right, 
the processes dominating the cooling curve are molec- 
ular hydrogen cooling (for T < 10*K), hydrogen line 
cooling (1st trough), helium line cooling (2nd trough), 
Bremsstrahlung from free electrons (for T > lO^K) and 
Compton cooling against cosmic microwave background 
photons (horizontal line for T > lO^K). The curve cor- 
responds to zero metallicity, since we are interested in 
whether the first galaxies can form. 

The cooling physics of course depends only on atomic 



^ This criterion is closely linked to the question of whether ob- 
servers can form if one is prepared to wait an arbitrarily long 
time. First of all, if the halo fails to contract substantially in a 
Hubble time, it is likely to lose its identity by being merged into a 
larger halo on that timescale unless x S> 1 so that p\ has frozen 
clustering growth. Second, if Tvir < 10'*K so that the gas is 
largely neutral, then the cooling timescale will typically be much 
longer than the timescale on which baryons evaporate from the 
halo, and once less than O.OSMq of baryons remain, star forma- 
tion is impossible. Specifically, for a typical halo profile, about 
1% of the baryons in the high tail of the Bolzmann velocity dis- 
tribution exceed the halo escape velocity, and the halo therefore 
loses this fraction of its mass each relaxation time (when this 
high tail is repopulated by collisions between baryons). In con- 
trast, the cooling timescale is linked to how often such collisions 
lead to photon emission. This question deserves more work to 
clarify whether all Tvir 3> lO^K halos would cool and form stars 
eventually (providing their protons do not have time to decay). 
However, as we will see below in Section IIVI this question is 
unimportant for the present paper's prime focus on axion dark 
matter, since once we marginalize over pA, it is rather the upper 
limit on density in Figure|21that affects our result. 
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processes, i.e., on the three parameters (a, /3, rrip). Re- 
quiring the cooling timescale to not exceed some fixed 
timescale would therefore give a curve independent of all 
cosmological parameters. Since we are instead requir- 
ing the cooling timescale ((x for all processes involv- 
ing particle-particle collisions) not to exceed the Hubble 

timescale (ex p-at^^ oc (n//b)~^/^), our cooling curve will 
depend also on the baryon fraction /;,, with all parts ex- 
cept the Compton piece to the right scaling vertically as 



2. Disk fragmentation and star formation 

We have now discussed how fundamental parameters 
determine whether dark matter halos form and whether 
gas in such halos can cool efficiently. If both of these 
conditions are met, the gas will radiate away kinetic en- 
ergy and settle into a rotationally supported disk, and 
the next question becomes whether this disk is stable or 
will fragment and form stars. As described below, this 
depends strongly on the baryon fraction /b = £,h/^, ob- 
served to be /b 1/6 in our universe. The constraints 
from this requirement propagate directly into Figure IT^ 
in Section llV Dl rather than Figure 13 

First of all, stars need to have a mass of at least M^i^ « 
O.OSMq for their core to be hot enough to allow fusion. 
Requiring the formation of at least one star in a halo of 
mass M therefore gives the constraint 



/b> 



M 



(36) 



An interesting point |lllj | is that this constraint places 
an upper bound on the dark matter density parame- 
ter. Since the horizon mass at equality is Moq ~ i^^^ 
|85| . the baryon mass within the horizon at equality is 
Meq£,h/^ ^ Cb/C^: which drops with increasing ^. This 
scale corr espo nds to the bend in the banana of Figure 13 
so unless |lllj | 



/b > M^U^ 



(37) 



one needs to wait until long after the first wave of halo 
formation to form the first halo containing enough gas 
to make a star, which in turn requires a correspondingly 
small pA-value. The ultimate conservative limit is re- 
quiring that ~ phH~^ , the baryon mass within the 
horizon, exceeds A/min. Since increases during mat- 
ter domination and decreases during vacuum domination, 
taking its maximum value ~ fhP^^^"^ when a; ~ 1, 
requirement gives 



the > M, 



fh Z MninPA ^■ 



(38) 



However, these upper limits on the dark matter param- 
eter density parameter S.c are very weak: for the ob- 
served values of ^b and p\ from Table 1, equations H37|) 



Cc < ^ ^h/M^inp]^"^ ~ 10~^, respectively, Hmits many 
orders of magnitude above the observed value ~ 10~^*. 

It seems likely, however, that Mmin is a gross underes- 
timate of the baryon requirement for star formation. The 
fragmentation instability condition for a baryonic disc is 
essentially that it should be self-gravitating in the "ver- 
tical" direction. This is equivalent to requiring the bary- 
onic density in the disc exceed that of the dark matter 
background it is immersed in. The first unstable mode is 
then the one that induces breakup into spheres of radius 
of order the disc thickness. This conservative criter ion i s 
weaker than the classic Toomre instability criterion |l02j| , 
which requires the disc to be self-gravitating in the radial 
direction and leads to spiral arm format ion. If the halo 
were a singular isothermal sphere, then |103| | instability 
would require 



/b> A 



1/2 



(39) 



Here Tmin is the minimum temperature that the gas 
can cool to and A is the dimensionless specific angu- 
lar momentum parameter, which has an approximately 
lognormal distribution centered around 0.08 ,.103] . For 
an upper limit on the dark matter parameter ^c, what 
matters is thus the upper limit on A, the upper limit 
on Tvir and the lower limit on Tmin, all three of which 
are quite firm. Ignoring probability distributions, taking 
Tmin a2/3TOp/61na-i - lO'^K - 1 eV (atomic Hydro- 
gen line cooling) and Tvir ~ 500eV (Milky Way) gives 
Tvir/Tmin 500, ^ SOO^b from equation Tak- 
ing the extreme values Tmin ~ 500K (_ff2-cooling freeze- 
out) and Tvir ^ m^Q ^ 20 keV (largest clusters) gives 
Tvh/Tniin ~ lOiog and ^ < lO^Qi/s^^^ _ iq*^^. Qu the 
other extreme, if we argue that Tvir ^ ?niin for the very 
first galaxies to form, then equation H39() gives ^/i^b ^ 12, 
which is interestingly close to what we observe. 

For NFW potentials 104], A-dependence is more com- 
plicated and the local velocity dispersion of the dark mat- 
ter near the center is lower than the mean Tvir. 

Even if the baryon fraction were below the threshold of 
equation H39|) , stars could eventually form because viscos- 
ity (even just "molecular" viscosity) would redistribute 
mass and angular momentum so that the gas becomes 
more centrally condensed. The condition then becomes 
that the mass of spherical blob of radius R ^ M^/T^i^ 
exceed the dark mass within that radius. For the 
isothermal sphere, this would automatically happen, but 
for a realistic flat-bottomed or NFW potential, then this 
gives a non-trivial inequality. For instance, in a parabolic 
potential well, the requirement would be that 



vir 



3/2 



(40) 



and ^ give < C ^ (Cb/Mmi„) 



1/3 



10 



-22 



and 



Although a weaker limit limit than equation (|39() , giving 
$c ^ lO'^^b for the above example with Tvir/Tmin ~ 500, 
it is still stronger than those of equation H37|) and equa- 
tion H38|l . A detailed treatment of these issues is beyond 
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the scope of the present paper; it should include mod- 
eling of the i — > CXI limit as well as merger-induced star 
formation and the effect of dark matter substructure dis- 
turbing and thickening the disk. 

C. Second generation star formation 

Suppose that all the above conditions have been met 
so that a halo has formed where gas has cooled and 
produced at least one star. The next question becomes 
whether the heavy elements produced by the death of the 
first star(s) can be recycled into a solar system around a 
second generation star, thereby allowing planets and per- 
haps observers made of elements other than hydrogen, 
helium and the trace amounts of deuterium and lithium 
left over from big bang nucleosynthesis. 

The first supernova explosion in the halo will release 
not only heavy metals, but also heat energy of order 

E = IQ-^rrip^ - lO^^J lO^^erg. (41) 

Here is the approximate binding energy of a Chan- 
drashekar mass at its Schwarzschild radius, with the pref- 
actor incorporating the fact that neutron stars are usu- 
ally somewhat larger and heavier and, most importantly, 
that about 99% of the binding energy is lost in the form 
of neutrinos. 

By the virial theorem, the gravitational binding energy 
E of the halo equals twice its total kinetic energy, i.e., 

E^Mvl,^^ J^S^, (42) 

V "l°Pvir 

where we have used equation (|23|l in the last step. If 
Esn » E, the very first supernova explosion will therefore 
expel essentially all the gas from the halo, precluding the 
formation of second generation stars. Combining equa- 
tions (|41|l and H42I) . we therefore obtain the constraint 

Pvir < 106^. (43) 
mp 

Figure|21illustrates this fact that lines of constant binding 
energy have slope 5, and shows that the second genera- 
tion constraint rules out an interesting part of the the 
(rivir, Tvir) plane that is allowed by both cooling and dis- 
ruption constraints. 

While we have ignored the important effect of cooling 
by gas in the supernova's immediate environment, this 
constraint is probably nonetheless rather conservative, 
and a more detailed calculation may well move it further 
to the right. First of all, many supernovae tend to go off 
in close succession in a star formation site, thereby jointly 
releasing more energy than indicated by equation (|42|l . 
Second, it is likely that many supernovae are required to 
produce sufficiently high metallicity. Since ~ 1 supernova 
forms per IOOMq of star formation, releasing ~ IM© 



of metals, raising the mean metallicity in the halo to 
solar levels 10~^) would require an energy input of 
order ^^/(lOOm^/mp) ~ lO^^mp ~ lOkeV per proton. 
A careful calculation of the corresponding temperature 
would need to model the gas cooling occurring between 
the successive supernova explosions. 

D. Encounters and extinctions 

The effect of halo density on solar system destruction 
was discussed in [85l | making the crude assumption that 
all halos of a given density had the same characteristic 
velocity dispersion. Let us now review this issue, which 
will play a key role in determining predictions for the 
axion density, from the slightly more refined perspective 
of Figure|21 Consider a habitable planet orbiting a star of 
mass M suffering a close encounter with another star of 
mass , approaching with a relative velocity f -j- and an 
impact parameter b. There is some "kill" cross section 
(T|(Af, Af|, Uj) = Tvb^ corresponding to encounters close 
enough to make this planet uninhabitable. There are 
several mechanisms through which this could happen: 

1. It could become gravitationally unbound from its 
parent star, thereby losing its key heat source. 

2. It could be kicked into a lethally eccentric orbit. 

3. The passing star could cause disastrous heating. 

4. The passing star could perturb an Oort cloud in the 
outer parts of the solar system, triggering a lethal 
comet impact. 

The probability of the planet remaining unscathed for a 
time t is then e~'''t*, where the destruction rate is 7| is 
n^,tT|W-f appropriately averaged over incident velocities W| 
and stellar masses M and M| . Assuming that 

1/2 

oc Uvir oc r^jj. and o-j is independent of Uvir and Tvir, 
contours of constant destruction rate in FigurelHare thus 
lines of slope —1/2. The question of which such contour is 
appropriate for our present discussion is highly uncertain, 
and deserving of future work that would lie beyond the 
scope of the present paper. Below we explore only a 
couple of crude estimates, based on direct and indirect 
impacts, respectively. 

1. Direct encounters 

Lightman |112| has shown that if the planetary surface 
temperature is to be compatible with life as we know it, 
the orbit around the central star should be fairly circular 
and have a radius of order 

rau - a-5^p 3/2/3-2 ^ iQiim, (44) 

roughly our terrestrial "astronomical unit" , precessing 
one radian in its orbit on a timescale 

torb a"l''^/2„^-5/4^-3 _ q ^^^^ _ (45) 
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An encounter with another star with impact parameter 
1" ^ ''au has the potential to throw the planet into a highly 
eccentric orbit or unbind it from its parent star.® For n^,, 
what matters here is not the typical stellar density in a 
halo, but the stellar density near other stars, including 
the baryon density enhancement due to disk formation 
and subsequent fragmentation. Let us write 



(47) 



where A^* ~ ~ 10^^ is the number of protons 

in a star and the dimensionless factor parametrizes 
our uncertainty about the extent to which stars concen- 
trate near other stars. Substituting characteristic values 
Uvir ~ IQr'/nv' for the Milky Way and '-^ (Ipc)"'^ 
for the solar neighborhood gives a concentration factor 



fi, ~ 10 . Using this value, 



and (T| — 



excludes the dark shaded region to the upper right in 
Figure 12 if we require to exceed imin ~ 10^ years 

^ a^ag "^^^Z?"^, the lifetime of a bright star Q. It is of 
course far from clear what is an appropriate evolutionary 
or geological timescale to use here, and there are many 
other uncertainties as well. It is probably an overesti- 
mate to take w-f v^-n- since we only care about relative 
velocities. On the other hand, our value of is an under- 
estimate since we have neglected gravitational focusing. 
The value we used for is arguably an underestimate as 
well, stellar densities being substantially higher in giant 
molecular clouds at the star formation epoch. 



2. Indirect encounters 

In the above-mentioned encounter scenarios 1-3, the 
incident directly damages the habitability of the planet. 
In scenario 4, the effect is only indirect, sending a hail 
of comets towards the inner solar system which may at 
a later time impact the planet. This has been argued to 
place potentially stronger upper limits on Q than direct 
encounters Isij . 

Although violent impacts are commonplace in our par- 
ticular solar system, large uncertainties remain in the 
statistical details thereof and in the effect that chang- 
ing and W| would have. It is widely believed that 



Encounters have a negligible effect on our orbit if they are adi- 
abatic, i.e., if the impact duration r/v 2> iorb so that the solar 
system returned to its unperturbed state once the encounter was 
over. Encounters are adiabatic for 



< 



1/4 



0.0001 ~ 30km/s, 



(46) 



the typical orbital speed of a terrestrial planet. As long as the 
impact parameter < rau, however, the encounter is guaranteed 
to be non-adiabatic and hence dangerous, since the infalling star 
will be gravitationally accelerated to at least this speed. 



solar systems are surrounded by a rather spherical cloud 
of comets composed of ejected leftovers. Our own par- 
ticular Oort cloud is estimated to contain of order 10^^ 
comets, extending out to about a lightyear 10^ AU) 
from the Sun. Recent estimates suggest that the impact 
rate of Oort cloud comets exceeding 1 km in diameter 
is between 5 and 700 per million years jll.ll |. These im- 
pacts are triggered by gravitational perturbations to the 
Oort cloud sending a small fraction of the comets into 
the inner solar system. About 90% of these perturba- 
tions are estimated to be caused by Galactic tidal forces 
(mainly related to the motion of the Sun with respect to 
the Galactic midplane), with random passing stars be- 
ing responsible for most of the remainder and random 
passing molecular clouds playing a relatively minor role 

jiia..ii4, .nli. 

It is well-known that Earth has suffered numerous vi- 
olent impacts with celestial bodies in the the past, and 
the 1994 impact of comet Shoemaker-Levy on Jupiter 
illustrated the effect of comet impacts. Although the 
nearly 10 km wide a steroid that hit the Yucatan 65 mil- 
lion years ago 'll(l| may actually have helped our own 
evolution by eliminating dinosaurs , a l arger impact of a 
30 km object 3.47 billion years ago |ll7l | may have caused 
a global tsunami and massive heating, killing essentially 
all life on Earth (For comparison, the Shoemaker-Levy 
fragments were less than 2 km in size.) We therefore 
cannot dismiss out of hand the possibility that we are 
in fact close to the edge in parameter space, with only a 
modest increase in comet impact rates on planets causing 
a significant drop in the fraction of planets evolving ob- 
servers. This possibility is indicated by the light-shaded 
excluded region in the upper right of Figure |21 

However, there are large uncertainties here of two 
types. First, we lack accurate risk statistics for our own 
particular solar system. Although the lunar crater radius 
distribution is roughly power law of slope —2 at high end, 
we still have very li mited knowledge of the size distribu- 
tion of comet nuclei |118| and hence cannot accurately es- 
timate the frequency of extremely massive impacts. Sec- 
ond, we lack accurate estimates of what would happen in 
denser galaxies. There, the more frequent close encoun- 
ters with other stars could rapidly strip stars of much 
of their dangerous Oort cloud, so it is far from obvious 
that the risk rises as n^t)f . One interesting possibility is 
that the inner Oort cloud at radii < lOOOAU contributes 
a substantial fraction of our impact risk, in which case 
such tidal stripping of most of the cloud by volume will 
do little to reduce risk. 

An indirect hint that comet impacts are anthropically 
important may be the observation |ll9j that the orbit of 
our Sun through the galaxy appears fine-tuned to min- 
imize Oort cloud perturbations and resulting comet im- 
pacts: compared to similar stars, its orbit has an un- 
usually low eccentricity and small amplitude of vertical 
motion relative to the Galactic disk. 
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3. Nearby explosions 

A final category of risk that deserves furtlicr explo- 
ration is that from nearby supernova explosions and 
gamma-ray bursts. Since these are independent of stel- 
lar motion and thus depend only on n*, not on Tvir, they 
would correspond to a horizontal upper cutoff in FigureEl 
For example, the Ordovician extinction 440 million years 
ago has been blamed on a nearby gamma-ray burst. This 
would have depleted the ozone layer causing a massive in- 
crease in ultraviolet solar radiation exposure and could 
also have triggered an ice age |12C| . 

E. Black hole formation 

There are two potentially rather extreme selection ef- 
fects involving black holes. 

First, Figure 121 illustrates a vertical constraint on the 
right side corresponding to large Tvir-values. Since the 
right edge of the halo banana is at Wvir ~ Q^^^c, typical 
halos will form black holes for Q- values of order unity as 
discussed in [s^. Specifically, typical fluctuations would 
be of black-hole magnitude already by the time they en- 
tered the horizon, converting some substantial fraction 
of the radiation energy into black holes shortly after the 
end of inflation and continually increasing this fraction 
as longer wavelength fluctuations entered and collapsed. 
For a scale-invariant spectrum, extremely rare fluctua- 
tions that are standard deviations out in the Gaus- 
sian tail can cause black hole domination if [s^ 

Q Z 10^^ (48) 

This constraint is illustrated in Figure E| below rather 
than in Figure |21 

A second potential hazard occurs if halos form with 
high enough density that collapsing gas can trap photons. 
This makes the effective 7-factor close to 4/3 so that the 
Jeans mass does not fall as collapse proceeds, and col- 
lapse proceeds in a qualitatively different way because 
there will be no tendency to fragment. This might lead 
to production of single black hole (instead of a myriad 
of stars) or other pathological objects, but what actu- 
ally happens would depend on unknown details, such as 
whether angular momentum can be transported outward 
so as to prevent the formation of a disk. 

F. Constraints related to the baryon density 

To conclude our discussion of selection effects, FigureEl 
illustrates a number of constraints that are independent 
of Q and pA, depending on the density parameters 
and for baryonic and dark matter either through their 
ratio or their sum ^ = + Cc- 

As we saw in Section llll B 21 a very low baryon fraction 
/b = ^h/S. ^ 300 may preclude disk fragmentation and 
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FIG. 6: Qualitatively different parts of the {^b, ?c)-plane. Our 
observed universe corresponds to the star, and altering the param- 
eters to leave the white region corresponds to qualitative changes. 
The hatched regions delimit baryon fractions fb = ^b/?c > 300 
and /fe < 1, which may suppress galaxy formation by impeding 
disk fragmentation and by Silk damping out Galaxy-scale mat- 
ter clustering, respectively, The shaded regions delimit interesting 
ranges of constant total matter density parameter 5 = 5b+?c — to 
what extent these qualitative transitions are detrimental to galaxy 
formation deserves further work. 



star formation. On the other hand, a baryon fraction 
of order unity corresponds to dramatically suppressed 
matter clustering on Galactic scales, since Silk damping 
around the recombination epoch suppresses fluctuations 
in the baryon component and the fluctuations that cre- 
ated the Milky Way were preserved through this epoch 
mainly by dark matter [90| . 

Recombination occurs at T ~Ry/50 ^ mpa^/3/100 ^ 
3000K, whereas matter-domination occurs at Tcq ^ 
0.2^. If <^ < 0.05mpQ;^/3, recombination would therefore 
precede matter-radiation equality, occurring during the 
radiation-dominated epoch. It is not obvious whether 
this would have detrimental effects on galaxy formation, 
but it is interesting to note that, as illustrated in FigureEl 
our universe is quite close to this boundary in parameter 
space. 

If we instead increase ^, a there are two qualitative 
transitions. 

In the limit x = pA/pm 1, using equation H15|l . 
equation l|A4|l and pm — S.rij gives 

" " (2^) ^ ^ 0.33028 s(m)^. (49) 



19 



This means that the first Galactic mass structures 
(s(/i) « 28) go nonUnear at T ~ Q^Q, i.e., before recom- 
bination if ^ ^ lO^^mpa^Pmp/Q. A baryonic cloud able 
to collapse before or shortly after recombination while 
the ionization fraction remained substantial would trap 
radiation and, 



as mentioned above in Section llll El po- 
tentially produce a single black hole instead of stars. 

Big Bang Nucleosynthesis (BBN) occurs at T ^ 
rripPn ^ 1 MeV, so if we further increase f so that 
^ rripPm then the universe would be matter-dominated 
bef ore B BN, producing dramatic (but not necessarily fa- 
tal |l05j |) changes in primordial element abundances. 



IV. PUTTING IT ALL TOGETHER 

As discussed in the introduction, theoretical predic- 
tions for physical parameters are confronted with obser- 
vation using equation where neither of the two factors 
/prioi(p) and /seiec(p) are optional. Above we discussed 
the two factors for the particular case study of cosmol- 
ogy and dark matter, covering the prior term in SectionllTI 
and the selection effect term in Section UTTl Let us now 
combine the two and investigate the implications for the 
parameters 5, Q and pA- 



A. An instructive approximation 

For this exercise, we would ideally want to compute, 
say, the function /sgIcc(p) defined as the fraction of pro- 
tons ending up in stable habitable planets as function 
of (pa, C; Cb, Q), leaving the particle physics parameters 
fixed at their observed values. Figure El shows that this 
is quite complicated, for two reasons. First, as discussed 
above we have computed only certain integrated versions 
of /scioc- Second, the other selection effects discussed 
involve substantial uncertainties. To provide useful qual- 
itative intuition, let us therefore start by working out 
the implications of the simple approximation that there 
are sharp upper and lower halo density cutoffs, i.e., that 
observers only form in halos within some density range 
Pmin < Pvir < Pmax, whcrc thcse two density limits may 
depend on Q, ^ and ^b- Based on the empirical obser- 
vation that typical galaxies lie near the right side of the 
cooling curve in Figure|21 one would expect Pmin to be de- 
termined by the bottom of the cooling curve and Pmax by 
the intersection of the cooling curve with the encounter 
constraint. 

Equation H2U|I showed that the fraction of 
all protons collapsing into a halo of mass fi is 
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Of 



where = Q i can 

be crudely interpreted as the characteristic density 
of the first halos to form. Roughly, A/s{^)Gao ^ 1, 
is the matter-radiation equality density and is 
the factor by which our Universe gets diluted between 
equality (when fluctuations start to grow) and galaxy 
formation. More generally, of these protons, the fraction 
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FIG. 7: The fraction of all protons that end up in halos with 
density in the range Pmin < Pvir < Pmax is shown as a function 
of the cosmological parameters pA and p» = The con- 

tours show fractions 0.3, 0.05 and 0.003 for the example Pmin = 



ru , Pmax 



10" 



(the two dashed lines) and ^ 



10" 



(corresponding roughly to IO^'^Mq halos and giving s(p) 28). 
When Pa S> Q^^*, essentially no halos at all are formed. When 
< lO^'^Pmini essentially all halos have Pvir < Pmin- When 
Q^?" > 102 

Pmax, essentially all halos have Pvir ^ Pmax- The star 
shows our observed parameter values. 



fh in halos within a density range pi 



/h(pA,C,Cb,(3) = erfc 



GA(a;(pmin)) 



-erfc 



(pa/p*) 



GA(a;(Pmax)) 

(50)" 

where the function x{p) is given by equation (|24|l ex- 
tended so that x{p) = oo for p < 16pA. For convenience, 
we have here defined the rescaled density 



A 



(51) 



Figure[7|is a contour plot of this function, and illustrates 
that in addition to classical constraint pa <^ Ql^ dat- 
ing back to Weinberg and others [HlM |m ITosl 

109, HQj, we now have two new constraints as well: pA- 
independent upper and lower bounds on p* = Q^^'^. 



B. Marginalizing over pA 

As mentioned in Section ^ there is good reason to 
expect the prior fpiior{pA,^,Q) to be independent of 
Pa across the tiny relevant range |pa| ^ Q^£,'^ where 
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FIG. 8: The probability factor fhip*) after marginalizing over p\. 
The black line of slope 1 is for Pniin = Oi Pmax = oo. Increas- 
ing Pmin shifts the exponential cutoff on the left side progressively 
further the right: p^in = 10"^, lO"", 10"^ and 10"^ times the 
density at the dotted horizontal line gives the four curves on the 
left. Decreasing pmax shifts left the break to a —1/3 slope further 
to the right: pmax = 100, 10, 1 and 0.1 times the density at the 
dotted horizontal line gives the four downward-sloping curves on 
the right, respectively. 



/scicc(/3A, C> Q) is non-negligible. This implies that we 
can write the predicted parameter probability distribu- 
tion from equation ^ as 

/(pa, Cb, Q) OC /prior(^, ^h,Q)fo{(„ Cb, Q)fh{PA, £.h, Q), 

(52) 

where fo{^,£,h,Q) is the product of all other selection 
effects that we have discussed — we saw that these were 
all independent of pA- 

The key point here is that the only pA-dependence in 
equation (|52|l comes from the /;i-term, i.e., from selec- 
tion effects involving the halo banana. It is therefore in- 
teresting to marginalize over pA and study the resulting 
predictions for Q and ^. Thus integrating equation ifS^ 
over PA, we obtain the theoretically predicted probability 
distribution 

fitQ) a /prior(C,'Cb,Q)/o(C,Cb,Q)//i(C,Cb,Q), where 



A(e,eb,Q) 



fhipA,t^h,Q)dpA. 



(53) 



Note that if pmin and pmax are constants, then Q) is 
really a function of only one variable: since equation H50|) 
shows that ^ and Q enter in fhip a, ^, Q) only in the com- 
bination p^, = Q"^^^, the marginalized result fh{£,,Q) 



will be a function of p* alone. Figure |S1 shows this 
function evaluated numerically for various values of Pmin 
and Pmax, as well as the following approximation (dotted 
curves) which becomes quite accurate for p <C Pmin and 

P > Pmax: 



724.9 



.1/3 
P* 

Pmax 



Go 



,erfc 



(tIP^) 



1/3- 



(54) 

This approximation is valid also when Pmin and Pmax are 
functions of the cosmological parameters, as long as they 
are independent of pA. To understand this result qual- 
itatively, consider first the simple case where we count 
baryons in all halos regardless of density, i.e., Pmin = 0, 
Pmax — oo. Then a straightforward change of variables 
in equation (|52l) gives fh{p*) oc p^, — Q'^S,'^, as was pre- 
viously derived in [a. This result troubled the authors 
of both and |8ll |. since any weak prior on Q from 
inflation would be readily overpowered by the Q'^-factor, 
leading to the incorrect prediction of a much larger CMB 
fluctuation amplitude than observered. The same result 
would also spell doom for the axion dark matter model 
discussed in Section Hi Al since the ^^^/^-prior would be 
overpowered by the ^^-factor from the selection effect and 
dramatically overpredict the dark matter abundance. 

Figure 13 shows that the presence of selection effects on 
halo density has the potential to solve this problem. Both 
the problem and its potential resolution can be intuitively 
understood from Figure[7| Note that fh{p*) is simply the 
horizontal integral of this two-dimensional distribution. 
For Pmin <C p* <C Pmax, the integrand ~ 1 out to the 
P* ^ PA-Iine marked "NO HALOS", giving fh{p*) oc p*. 
In the limit p* 3> Pmax, however, we are integrating 
across the region marked 'TOO DENSE" in Figure|Hl and 

the result drops as fh{p*) oc pmax/pi^^ = Pmax/QC''^^- 
This is why Figure El shows a break in slope from -1-1 
to —1/3 at a location p^, ~ pmax/200, as illustrated by 
the four curves with different pmax- values. Conversely, in 
the limit p^ <^ Pmin, we are integrating across the expo- 
nentially suppressed region marked 'TOO DIFFUSE" in 
Figure causing a corresponding exponential suppres- 
sion in Figure 13 for p» ^ Pmin- 

Figures 1^ and [TUI illustrate this for the simple example 
of using the Section FlI Al axion prior /prior(0 
treating Q as fixed at its observed value, ignoring addi- 
tional selection effects (/o(C) = 1), and imposing a halo 
density cutoff Pmax of 5000 times the current cosmic mat- 
ter density. 

The key features of this plot follow directly from our 
analytic approximation of equation (|54|l . First, for p, ^ 
Pmax, we saw that //i oc p* oc so equation H53II gives 
the probability growing as / cx ^^^Z^^'*!^ — <^^'^ (the 
last (^-factor comes from the fact that we are plotting 
the probability distribution for log^ rather than ^). This 
means that imposing a lower density cutoff Pmin would 
have essentially no effect. This also justifies our approxi- 
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FIG. 9: Probability distribution for the axion dark matter den- 
sity parameter J measured from a random 10^^ Mq lialo with 
virial density below 5000 times the present cosmic matter den- 
sity. Green/light shading indicates the 95% confidence interval. 
The dotted vertical line shows our observed value ^ ss 4 eV, in 
good agreement with the prediction. The dashed curves show the 
analytic asymptotics ^^/^ and respectively. 



mation of using the total matter density parameter ^ as a 
proxy for the dark matter density parameter S^^- we have 
a baryon fraction ^b/fc <C 1 in the interesting regime, 
with a neghgible probabihty for Cc ^ ~ 0.6eV (the 
probabihty curve continues to drop still further to the 
left even without the /^.-factor). Second, for ^ Pmax, 
we saw that fh oc p7^^^ esc so equation l(5^ gives 

the probability falling as /(In^ cx ^^^^^^^^^ = 

Although this simple example was helpful for build- 
ing intuition, a more accurate treatment is required be- 
fore definitive conclusions about the viability of the axion 
dark matter model can be drawn. One important issue, 
to which we return below in Section Hvl is the effect of 
the unknown Q-prior, since the preceding equations only 
constrained a combination of Q and ^. A second impor- 
tant issue, to which we devote the remainder of this sub- 
section, has to do with properly incorporating the con- 
straints from Figure 13 If we only consider the encounter 
constraint and make the unphysical simplification of ig- 
noring the effect of halo velocity dispersion, then our up- 
per limit should be not on dark matter density (pmax 
constant) but on baryon density (pmax'^b/C constant). In- 
serting this into equation (|54|l makes the term p^'^/pfnax 
independent of ^ as ^ ^ oo, replacing the fh oc Q^^^^**/^ 

^ 4/3 

cutoff by ffi oc Q~^y^ , which is constant if ^b is. The 
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FIG. 10: Same as previous figure, but showing the probability 
distribution for the quantity which has a uniform prior in our 

axion model. The shape of this curve therefore reflects the selection 
cfi'ect alone. 



axion prior would then predict a curve /(In^) oc i^^/^ ris- 
ing without bound. This failure, however, results from 
ignoring some of the physics from Figure 13 Consider- 
ing a series of nominally more likely domains with pro- 
gressively larger ^, they typically have more dark energy 
(pa oc ^) and higher characteristic halo baryon density 
(nvir oc ^^£,h), so stable solar systems are found only in 
rare galaxies that formed exceptionally late, just before 
PA-domination, with rtvir oc Pvir^b/C ^ ^^Ph^h/^ oc ,^b 
roughly ^-independent and below the the maximum al- 
lowed value. Since the increased dark matter density 
boosts virial velocities, this failure mode corresponds 
to moving from the star in Figure |21 straight to the 
right, running right into the constraints from cooling and 
velocity-dependent encounters. In other words, a more 
careful calculation would be expected to give a prediction 
qualitatively similar to that of Figure^ The lower cutoff 
would remain visually identical to that of FigureElas long 
as Pmin can be neglected in the calculation, whereas the 
upper cutoff would become either steeper or shallower 
depending on the details of the encounter and cooling 
constraints. This interesting issue merits further work 
going well beyond the scope of the present paper, ex- 
tending our treatment of halo formation, halo mergers 
and galaxy formation into a quantitative probability dis- 
tribution in the plane of Figure |2| i.e., into a function 
/(Tvir, nviri Q, Pa, ?b) that could be multiplied by the 
plotted constraints and marginalized over (Tvir,rivir) to 
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FIG. 11: Probability distribution for the quantity R = pa/^'^Q^ 
measured from a random 10^^ Mq halo, using a uniform prior for R 
and ignoring other selection effects. This is equivalent to treating ^ 
and Q as fixed. Green/light shading indicates the 95% confidence 
interval, the dotted line indicates the observed value R 15. 



the theory predicts a 2-dimensional distribution in the 
[PA, 0~P^^^^ of which Figure 1^1 is the marginal distribu- 
tion for ^. The corresponding marginal distribution for 
Pa (marginalized over ^) will generally differ from Fig- 
ure llll in both its shape and in the location of its peak. 
It will differ in shape because the (pA, C)"distribution is 
not generally separable: the selection effects (as in Fig- 
ure lllfl will not be a function of pA times a function of ^; 
in other words, a uniform prior on pA will not correspond 
to a uniform prior on the quantity R = pa/Q^^^ plotted 
in Figure [7| once non-separable selection effects such as 
ones on the halo density parameter p* = Q^^* are in- 
cluded. Moreover, the distribution will in general not 
peak where that in Figure [TTI does because the predicted 
magnitude pA ''^ 10^^^'^ no longer comes from condition- 
ing on astrophysical measurements of Q and ^, but from 
other parameters like a, f3 and nip that determine the se- 
lection effects in Figure Inland the maximum halo density 
Pmax- In this case, whether the observed pA- value agrees 
with predictions or not thus depends sensitively on how 
strong the selection effects against dense halos are. 

In summary, the prediction of pA given Q and ^ is an 
unequivocal success, whereas predictions for Q, ^ and the 
entire joint distribution for (pa, ^, Q) are fraught with the 
above-mentioned uncertainties. Since a uniform pA-prior 
does not imply a uniform R = pA/Q'^^^-prior, we cannot 
conclude that anthropic arguments succeed in predicting 
R = pa/Q^£,^ [13 without additional hypotheses. 



give fscicciQ,PA,£.,^h)- 



C. Predictions for pA 

This above results also have important implications 
for the dark energy density pA. The most common 
way to predict its probability distribution in the liter- 
ature has been to treat all other parameters as fixed 
and assume both that the number of observers is propor- 
tional to the matter fraction in halos (which in our no- 
tation means /scioc(p) = erfc[A(pA/Q^$^)^/^/s(p)G'oo]) 
and that /prior (pa) is constant across the narrow range 
IpaI ^ Q^^'^ where /scioc is non-negligible. This gives the 
familiar result shown in Figure lTTI a probability distribu- 
tion consistent with the observed pA-value but favoring 
slightly larger values. The numerical origin of the pre- 
dicted magnitude pA ~ 10~^^^ is thus the measured value 

of p, = g3^4 

Generalizing this to the case where Q and/or ^ can vary 
across an inflationary multiverse, the predictions will de- 
pend on the precise question asked. Figure lTTl then shows 
the successful prediction for pA given (conditionalized on) 
our measured values for Q and ^. However, when test- 
ing a theory, we wish to use all opportunities that we 
have to potentially falsify it, and each predicted param- 
eter offers one such opportunity. For our axion example. 



D. Constraint summary 

Table 4 summarizes the constraints that we have dis- 
cussed above. Those not superseded by other stronger 
ones are also illustrated in Figure ^1 with i^b fixed at its 
observed value. Let us briefly comment on how they all 
fit together. 

The first key point to note in Figure ^1 is that the 
various constraints involve many different combinations 
of Q and ^, thereby breaking each other's degeneracies 
and providing interesting constraints on both ^ and Q 
separately. The fragmentation constraint (requiring disk 
instability) constrains ^ alone when is given whereas 
there are constraints on Q alone from both line cool- 
ing freezeout and primordial black hole trouble. The en- 
counter and cooling constraints both have different slopes 
(— 1 and —2/3, respectively) because they involve differ- 
ent powers of the baryon fraction. Whereas the non- 
linear halo constraint depends on the total halo density 
oc Q^S,*, the encounter constraint depends on the halo 
baryon density oc Qff^^b- The cooling constraint (which 
is equation (11) in [83 after changing variables to reflect 
the notation in this paper) comes from equating the cool- 
ing timescale (given by the inverse of the baryon density 
oc Q^£,'^£,h) and the dynamical timescale (given by the to- 
tal density as oc p^ir^^ c>c (Q'^^'*)^^^^, thus constraining 
the combination Q^£,^£,^- 

The second key point in Figure IT^ is that this combi- 
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Table 4: These constraints (see text) are summarized in Figure IT2I 



Constraint 


Generally 


Fixing (a, /3, r?ip) 


Fixing all but (Q, 


Need nonlinear halos 


IpaI £ P* 


IpaI/?*Q^ ^ 1 


Q > io-5(«/5o)-*/3 


Avoid line cooling freezeout 




Q > 10-** 


Q > 10-8 


Primordial black hole excess 


Q < 10-1 


Q < 10-1 


Q < 10-1 


Need cooling in Hubble time 


Q'i^^e > a-3ln[a-2]-l6/3/34m6/125 


Q3^2g2 > 10-129 


Q > io-8(5/Co)-2/3 


Avoid close encounters 




Q^'ib^ < 10-12=* 


Q < io-*(«/Co)-i 


Go nonlinear after decoupling 




< 10-30 


Q < io-2(€/Co)-i 


Need equality before decoupling? 


^ > 0.05a2/3mp 


« > 10-28 


(5/«o) > 1/3 


Avoid severe Silk damping 


/b < 1/2 


«/Cb > 2 


mo) > 1/3 


Need disk instability 


/b < 102 


«c/Cb < 102 


(S/Co) < 20 



cy io-» f 



Black hole trouble 




^, Matter density per CMB photon [eV] 

FIG. 12: Crude summary of constraints Q)-plane constraints 
from Table 4. The star shows our observation Q) ft: (4eV, 2 X 
IQ- ^). The parallel dotted lines are lines of constant characteristic 
halo density, so the constraint that dark matter halos form (that 
they go nonlinear before vacuum domination freezes fluctuation 
growth) rules out everything to the lower left of these lines, which 
from left to right correspond to Pa/Pa'"^ = 10-3, 1, 103. This 
means that if § is fixed, marginalizing over pA pushes things to the 
triangle at larger Q (the "smoothness problem" 1^), but if ^ too 
can vary, marginalizing over pA instead pushes things to the square 
at smaller Q. 



nation of constraints reverses previously published con- 
clusions about Q. Where we should expect to find our- 
selves in Figure El depends on the priors for both Q and 
Recall that halos form only above whichever dotted 
line corresponds to the pA-value. If the prior for pA is 
indeed uniform, then it is a priori (before selection ef- 
fects are taken into account) 1000 times more likely for 
the halo constraint to be the upper dotted line than the 
middle one, and this is in turn 1000 more likely than 
the bottom dotted line. In other words, marginalizing 



over Pa relentlessly pushes us towards the upper right 
in Figure ^1 towards larger value of p* = Q^£,'^ ■ This 
means that if ^ is fixed and only Q can vary across the 
ensemble, marginalizing over pA pushes things towards 
the triangle at larger Q. This "smoothness problem " 
was pointed out in 85] and elaborated in [gI I121L [l22j| . 
suggesting that unless the /prior (Q) falls off as Q"'* or 
faster, it is overpowered by the selection effect and one 
predicts a clumpier universe than observed (the star in 
Figure El should be shifted vertically up against the edge 
of the encounter constraint, to the triangle). Figure IT^ 
shows that if ^ too can vary, marginalizing over pA in- 
stead pushes things towards smaller Q, to the square, 
since this gives the greatest allowed halo density oc Q^£,'^- 
In , it was found that certain low-energy inflation sce- 
narios naturally predicted a rather flat priors for Ig Q — 
to determine whether they are ruled out or not thus re- 
quires more detailed modeling of the effect of the dark 
matter density on galaxy formation and encounters. 

Table 4 illustrates that (except for one rather unim- 
portant logarithm), all constraints we have considered 
correspond to hyperplanes in the 7-dimensional space 
parametrized by 



(Ig a, Ig P, Ig mp, Ig PA, Ig Q, Igt Ig Cb), 



(55) 



providing a simple geometric interpretation of the selec- 
tion effects. This means that if all priors are power laws, 
the preferred parameter values are determined by solving 
a simple linear programming problem, and will generi- 
cally correspond to one of the corners of the convex 7- 
dimensional allowed region. Above we have frequently 
used such sharp inequalities to help build intuition for the 
underlying physics. For future extensions of this work, 
however, it should be borne in mind that such inequalities 
are not necessarily sufficient to capture the key features 
of the problem, and that they can either overstate or 
understate the importance of selection effects. The ulti- 
mate goal is to test theories by computing the predicted 
probability distribution from equation , and almost no 
selection effects /soicc(p) correspond to sharp cutoffs re- 
sembling a Heaviside step function. Some (like the halo 
constraint) are exponential and thus able to overpower 
any power law prior. Others, however, may be softer, and 
it is therefore crucial to check whether their functional 
form falls of faster than than the relevant prior grows — 
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if not, the "constraint" will be penetrated, and the most 
likely parameter values may lie in the allegedly "disal- 
lowed region". Conversely, if /(p) = /prior(p)/scicc(p) 
is rather flat in a large "allowed region", then the most 
likely parameter values can lie well inside it, far from any 
edges, since it corresponds to the peak of the probability 
distribution rather than its edge. 

V. CONCLUSIONS 

We have discussed predictions for the 31-dimensional 
vector p of dimensionless physical constants in the con- 
text of "ensemble theories" producing multiple Hubble 
volumes ( "universes" ) between which one or more of the 
parameters can differ. The prediction is the probabil- 
ity distribution of equation |^ , where neither of the two 
factors is optional. Although it is generally difficult to 
evaluate either one of the terms, it is nonetheless cru- 
cial: if candidate theories involve such ensembles, they 
are neither testable nor falsifiable unless their probabilis- 
tic predictions can be computed. 

Although many scientists hope that /prior in equa- 
tion will be a multidimensional Dirac (5-function, ren- 
dering selection effects irrelevant, to elevate this hope 
into an assumption would, ironically, be to push the 
anthropic principle to a hedonistic extreme, suggesting 
that nature must be devised so as to make mathematical 
physicists happy. We must therefore face some difficult 
questions both about what, in principle, to select on, 
and how to quantitatively calculate that selection effect. 
For both statistical mechanics and quantum mechanics, 
it proved challenging to correctly predict probability dis- 
tributions. If the analogous challenge can be overcome 
for theories predicting dimensionless constants, they too 
will become testable in the same fashion. The fact that 
we can only observe once is not a show-stopper: a sin- 
gle observation of a vertically polarized photon passing 
through a polarizer 6 = 89.9° from vertical would rule out 
quantum mechanics at 1 — cos^ 9 k, 99.9997% confidence. 
The outstanding difficulty is therefore not a "philosoph- 
ical" problem that we cannot observe the full ensemble, 
but rather that we generally do not know how to compute 
the theoretical prediction /(p). 



A. Axion dark matter 

We have tackled this problem quantitatively for a par- 
ticular example where both factors in equation Q o,re 
computable: that of axion dark matter with its phase 
transition well before the end of inflation. We found 
that /prior(0 ^ that the prior dark mat- 

ter distribution has no free parameters even though the 
axion model itself does. We saw that this useful pre- 
dictability gain has the same basic origin as that of 
the zero-parameter (flat) pA-prior of [s^: the anthrop- 
ically relevant range over which /sciec is non-negligible 



is much smaller than the natural range of /prior, so 
only a particular property of /prior matters. For the fa- 
mous pA-case, the natural /prior-range involves either the 
supersymmetry-breaking scale or the Planck scale |pa| ^ 

Pplanck 

whereas the relevant range is |pa| ^ 10 ^^^Ppianck, 
so the property that /prior (pa) is smooth on that tiny 
scale implies that it is for all practical purposes constant. 
For our ^-case, the natural /prior-range is < ^ < ^* 

where lies somewhere between the inflation scale 
and the Planck scale, whereas the relevant range is much 
smaller, so the only property of /prior (^) that matters is 
its asymptotic scaling as ^ — > 0. 

In computing the selection effect factor /scioc(pAj Q)i 
we took advantage of the fact that none of these pa- 
rameters affect chemistry or biology directly. We could 
therefore avoid poorly understood biochemical issues re- 
lated to life and consciousness, and could limit our cal- 
culations to astrophysical selection effects involving halo 
formation, galaxy, solar system stability, etc. The same 
simplifying argument has been previously applied to pA 
and the neutrino density {e.g., [S^). We discovered that 
combining various astrophysical selection effects with all 
three parameters (pA^QA) ^s variables reversed previ- 
ous conclusions fFigure lT^ , pushing towards lower rather 
than higher Q- values. In particular, we found that one 
can reach misleading conclusions by simply taking /sciec 
to be the matter fraction collapsing into halos, neglecting 
the fact that very dense halos limit planetary stability. In 
particular, we found that if we impose a stiff upper limit 
on the density of habitable halos, the predicted proba- 
bility distribution for the dark matter density parameter 
(Figure inj) can be brought into agreement with the mea- 
sured ^- value for reasonable assumptions, i.e., that the 
pre-inflationary axion model is a viable dark matter the- 
ory. 



B. Multiple dark matter components 

Suppose that supersymmetry were to be discovered at 
the Large Hadron Collider, indicating the existence of a 
weakly interacting massive particle (WIMP) with a relic 
density approximately equal to the observed dark matter 
density, i.e., with ^wimp ~ Cc- There would then be a 
strong temptation to declare the dark matter problem 
solved. Based on our results, however, that would be 
quite premature! 

As detailed in Section Hi Ij| it is not implausible that 
the function /prior (^wimp) to be taken as input to the as- 
trophysics calculation is a sharply peaked function whose 
relative width is much narrower than that of /scicc(^c), 
perhaps only a few percent. This situation could arise if 
the fundamental theory predicted a broad prior only for 
the ^^-parameter in Table 1 (which controls the Higgs 
vacuum expectation value), since this may arguably be 
determined to within a few percent from selection effects 
in the nuclear physics sector, notably stellar carbon and 
oxygen production. If the pre-inflationary axion model 
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is correct, we would then have 



wimp 



« Cxlon and C^, 



(56) 



would be for all 



where /prior 

practical purposes a constant. This gives /prior (Cc) oc 
(Cc - ^wimp)"^/^ for 5c > ^wirnp, zcro othcrwisc, and the 
testable prediction becomes this times the astrophysical 
factor selection effects factor /seiec(Cc) that we have dis- 
cussed. Repeating the calculation of this implies 
that we should expect to find roughly comparable densi- 
ties of WIMP and axion dark matter. 

As shown in 14], removing the assumption that 
/prior (Cwimp) is sharply peaked does not alter this qual- 
itative conclusion. This conclusion would also hold for 
any other axion-like fields that were energetically irrele- 
vant during inflation, even unrelated to the strong CP- 
problem. More generally, string-inspired model building 
suggests the possibility that multiple species of "GIMPS" 
(gravitationally interacting massive particles that do not 
couple through the electromagnetic, weak or strong in- 
teractions) may be present. According to recent under- 
standing, their presence may be expecte d for modular 
inflation, though not for brane inflation [l23j . If such 
particles exist and have generic not-too-extreme priors 
(falling no faster than and rising no faster than the 
selection effect cutoff), then our anthropic constraints on 
the sum of their densities again predict a comparable den- 
sity for all of species of dark matter that have a rising 
prior .14] . 

Even before the LHC turns on, more careful calcu- 
lation of the astrophysical selection effects determining 
/seicc could give interesting hints. If the constraints on 
dense halos turn out to be rather weak so that the im- 
proved version of Figure El predicts substantially more 
axion dark matter than observed, this will favor an alter- 
native model such as WIMPs where /prior (Cc) is so nar- 
row that is effectively determined non-anthropically. 
Alternatively, if Figure |21 remains consistent with obser- 
vation, then it will be crucial to follow up an initial LHC 
SUSY or other WIMP detection, if it occurs, with de- 
tailed measurements of the relevant parameters to estab- 
lish whether the WIMP density accounted for (which a 
linear collider might, in favorable cases, pinpoint within 
a few percent 73]) is exactly equal to the cosmologically 
measured dark matter density, or whether a substantial 
fraction of the dark matter still remains to be accounted 
for. 
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APPENDIX A: FITTING FUNCTIONS USED 

In this Appendix, we derive a variety of fitting func- 
tions used in the paper, explicitly highlighting the de- 
pendence of standard cosmological structure formation 
on the dimensionless parameters ^, pA, Q and Since 
some of our approximations are quite accurate, we retain 
the exact numerical coefficients where appropriate. 



1. The fluctuation time dependence G{x) 



As shown in 



the function 



Ga(x)«x1/3 



1 



where a = 159/200 = 0.795 and 



-l/3o 



Go 



5r(|)r(|) 



30F 



1.43728, 



(Al) 



(A2) 



describes how, in the absence of massive neutrinos, fluc- 
tuations in a Pa > ACDM universe grow as the cos- 
mic scale factor a as long as dark energy is negligible 
(GA(a:) « x^/^ oc a (X (1 -I- z)~^ for a; <C 1) and then 
asymptote to a constant value as t oo and dark en- 
ergy dominates (GA(a;) —t Goo as a; ^ oo). 

Including the effect of radiation domination early on 
but ignoring baryon-photon coupling, the total growth 
factor is well approximated by |83| ] 



G{x) « 1 + 



-x-^'^Gk{xI 



(A3) 



which is accurate to better than 1.5% for all x In 
essence, fluctuations grow as (5 oc a oc x^^^ between mat- 
ter domination (x = x^q) and dark energy domination 

{x = 1), giving a net growth of Xc^^^ ■ Equation (|A3(I 
shows that they grow by an extra factor of 1.5 by start- 
ing slightly before matter domination and by an extra 
factor Ga(oo) k, 1.44 by continuing to grow slightly after 
dark energy domination. 

Let us now simplify this result further. From the defi- 
nition of ^, the matter density can be written prn = 
where is the number density of photons. The latter 
is given (in Planck units) by the standard black body 
formula 



2C(3) 



T^ C(3) ~ 1.20206. 



(A4) 



26 



The corresponding standard formula for the energy den- 
sity is 



amphtude Sh on the horizon. As our measure of comov- 
ing scale, we use the dimensionless quantity 



for photons and 



Pi 



21 / 4 \ 



15^ 



rp4 



4/3 



p2 = -\-) P-<~ 0.68132P, 



8 Vll 



(A5) 



(A6) 



for three standard species of relativistic massless neutri- 
nos. The matter-radiation equality temperature where 
Pm ~ Pi + pZ is thus 



T — 

cq — 



30C(3) 



21/4 

^ + ¥11 



4/3- 



0.220189^, (A7) 



measured to be about 9400K in our Hubble volume. 
Since Xcq = pA/Pm — Ph/^n^{Tcq), we thus obtain 



^14 



2 • ZQ^CiiY 



4/3- 



8 Vll 



^« 384.554^ 

^4 ^4 



(A8) 

The above relations are summarized in Tables 2 and 3. 
In all cases of relevance to the present paper, the growth 
factor G{x) 3> 1, allowing initial fluctuations Q <C 1 to 
grow enough to go nonlinear. We can therefore drop the 
first term in equation (jA3|) . obtaining the useful result 



G{x) 



45 • 21/3(^(3)4/3 



.14/3 



21 /4\'^' 



f4/3 

0.206271^Ga(x). 
Pa 



^4/3 

Pa 

(A9) 



2. The fluctuation scale dependence s(/i) 

The r.m.s. fluctuation amplitude cr in a sphere of radius 
R are given by (106.] 



An 



a;3/3 



P(fc) 



(2^' 



(AlO) 



where x = kR. We are interested in the linear regime 
long after the relevant fluctuation modes have entered the 
horizon, when all modes grow at the same rate (assum- 
ing a cosmologically negligible contribution from massive 
neutrinos as per |89l | . This means that a can be factored 
as a product of a function of time x and a function of 
comoving scale: 



cr = QG{x)s{^). 



(All) 



Here Q is the amplitude of primordial fluctuations cre- 
ated by, e.g., inflation, defined as the scalar fluctuation 



M = CM, 



(A12) 



where M is the comoving mass enclosed by the above- 
mentioned sphere. 

The horizon mass at matter-radiation equality is of 
order i^^^, so /i can be interpreted as the mass rela- 
tive to this scale. This equality horizon scale is the 
key physical scale in the problem, and the one on which 
the matter power spectrum has the well-known break in 
its slope, (fluctuation modes on smaller scales entered 
the horizon before equality when they could not grow). 
For reference, the measured vales from WMAP-I-SDSS 
give ^ R:! 3 X 10~^^, « 10^^ solar masses and 

fi « R/ {60h~^Mpc)^ for a sphere of comoving radius R. 
Numerically, using equation IjAlOp and employing j93 to 
compute P{k), we find that the scale-dependence s(/x) is 
approximately given by 



s(m) « {9.1^i-'^^Y + (50.5 lg(834 + ^"^Z-^) - 92)'^ 

(A13) 

where (3 — —0.27. This assumes that the primordial 
fluctuations are approximately scale-invariant and that 
massive neutrinos have no major effect, as indicated by 
recent measurements [s^ |^ For our observed ^- 

value a galactic mass scale M = IO^^Mq thus corre- 
sponds to fi — ^"^M K, 10^^ and s(/i) « 28. Again using 
measured parameter values, the common reference scale 
R = 8/i~iMpc corresponds io ji ~ 0.002 and s(/u) « 11. 



-l/3x 



1//3 



3. The virial density pvir(x) 

We will now derive the virial density of a halo formed 
at time a; in a ACDM universe, generalizing the standard 
result pvij_~ 187r^pm(a;). Whereas past work on this 
subject [M. I95I I has focused on accurate numerical fits 
to the observationally relevant present epoch and recent 
past (x < 1), we need a formula accurate for all times x. 

Consider first the evolution of a top hat overdensity 
in a ACDM universe. By Birkhoff 's theorem, it is given 
by the Friedman equation corresponding to a closed uni- 
verse: 

37j2 ^ K 1 

^ = PA + P.n - ^ = (1 - ^ + ^)PA, (A14) 

where we have defined a dimensionless scale factor A = 
(pa/Piii)^^"^ and a dimensionless curvature parameter k = 
{A/a^k. One readily finds that this will give H — and 
recoUapse for any curvature k > Kmin = 3/2^/'^ ~ 1.8899, 
so the scale factor ^max(K) at turnaround is given by 
solving equation ljA14p for H — 0, i.e., by numerically 
finding the smallest positive root A of the cubic polyno- 
mial A^-kA+1 = 0. ^max(Kmin) = 2^1/^ ^ 0.79, falling 
off as ^max(K) w 1/k for 1. 
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The result is shown in Figure [TSl together with three an- 
alytic fitting functions. The simple fit 



Pturn(a;) 



97r2 



PA 



(A18) 



is accurate to better than 18% for all x and becomes 
exact in the two limits a; — > and a; — > oo. The fit 



Pturn(a;) 



9^ 

4x 



6.6 

^,0.3 



PA 



(A19) 



is accurate to better than 4% for all x and also becomes 
exact in both limits, but we will nontheless use the less 
accurate equation (|A18|I because it is analytically invert- 
ible and still accurate enough for our purposes. Apply- 
ing the virial theorem gives a characteristic virial density 
Oyir '-^ 8/9turn, with a Weak additional dependence on p\ 
This means that a halo that virialized at time x 

has 



{ISir^p^ix))^ + {16pA 



200 
107 107 
I 200 



(A20) 



FIG. 13: The top panel shows the turnaround (minimum) matter 
density pturn of a top hat halo as a function of the time x when 
its collapses. The virial density of the resulting dark matter halo 
is Pvir ~ Spturn- The Solid curve is the exact result, with the 
asymptotic behavior ptum — > 97r^pni/4a:: as a; — > and pturn —> 2pA 
as a; — > oo (dotted lines). The three analytic fits correspond to 
eq uati on lAlSi (short-dashed), equation IA19i (dotted) and that 
of 0^5], which in our notation is (IStt^/x -|- IStt^ - 82 - 39x/(l-|- 
x))/8. The first two are seen to be accurate to better than 18% 
and 4% respectively (bottom panel). 



The age of the universe tturn when the top hat turns 
around is given by integrating dt = din A /{din A/ dt) = 
dA/{AH) using equation (|A14|I . i.e., 



dA 



A^-n^\ 



(A15) 



where = (87rG'/9A/3)^/^. For the unperturbed back- 
ground ACDM universe where we used x as a more con- 
venient time variable, we have k = and A = x 
when the overtensity turns around, giving tturn in terms 

of Xturn: 



1/3 
turn 



u 



1/3 



dA 



(A16) 



Eliminating fturn between equations (|A15|I and HA16|I al- 
lows us to numerically compute the function ^(xturn)- 

Putting everything together now allows us to numeri- 
cally compute the density of our top hat overdensity at 
turnaround: Since = p^/A^, we have 



Pturn U 



PA 



(A17) 



i.e., essentially the larger of the two terms l^n"^ pni{x) 
and 16pA- Inverting this, we obtain 



( Pv 



- 1 



(A21) 



A familiar quantity in the literature is the collapse 
overdensity Ac = Pvir/Pm — Pvirx/ p^, equaling 187r^ « 
178 for a flat p^ ~ universe. Using equation ljA18|l . 
equation ljA19|) and the fit in to the calculations of 
gives 



Ac(x) 



187r^ 



187r^ + 52.8a;' 



8a; 
9^ 

,0.7 



107 
200 



2QQ 
107 



16a:, 



187r^ -f (187r^ - 82)a; - 39 



1 



(A22) 

(A23) 

(A24) 



respectively. Figure 1131 shows that our approximations 
substantially improve the accuracy, equation ljA23|l being 
good to 4%. This improvement not surprising, since the 
approximation (|A24ll was not designed to work well for 
a; > 1. 



^ The detailed calculation of 19381 shows that the standard collapse 
factor of 2 varies at thr 20% level with pA at the 20% level. 
Qualitatively, positive (repulsive) pA increases the collapse factor 
because shells have to fall further to acquire a velocity which will 
bring the system into equilibrium. We will not explicitly model 
this correction here given the larger uncertainties pertaining to 
other aspects of our treatment. 
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